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Abstract 


A  new  high  pressure  combustion  facility  was  built  to  investigate  mixing  in 
axisymmetric,  turbulent  jets  exiting  into  quiescent  reservoirs.  The  facility  uses 
fluorine  and  nitric  oxide,  diluted  with  nitrogen,  for  chemical  product  formation 
that  is  accompanied  by  heat  release.  The  average  temperature  was  measured  by  a 
set  of  long,  thin,  resistance  wire  thermometers  stretched  across  the  jet  centerline  at 
16  downstream  locations  from  x/da  =  30  to  240.  Runs  at  several  stoichiometric 
mixture  ratios  <j>,  for  Reynolds  numbers  ranging  from  10,000  <  Rt  <  150,000, 
were  performed  to  determine  any  dependence  of  flame  length  on  Reynolds  number. 
The  Reynolds  number  was  varied  through  density,  z.e.,  pressure,  while  the  jet  exit 
velocity  and  exit  diameter  were  held  constant.  The  time-averaged  line  integral  of 
temperature,  measured  along  the  transverse  axis  of  the  jet  by  the  wires,  displays 
a  logarithmic  dependence  on  r/d*  within  the  flame  zone,  and  asymptotes  to  a 
constant  value  beyond  the  flame  tip,  as  predicted  from  sczding  and  similarity  argu¬ 
ments  for  a  momentum-dominated,  turbulent  jet.  The  main  result  of  the  work  is 
that  the  flame  length,  zis  estimated  from  the  temperature  measurements,  varies  with 
changes  in  Reynolds  number,  suggesting  that  the  mixing  process  is  not  Reynolds 
number  independent  up  to  i?e  =  150,000.  Specifically,  the  normalized  flame  length 
Li/d*  displays  a  lineeir  dependence  on  <f>,  with  a  slope  that  decreases  from  Re  = 
10,000  to  20,000,  and  then  remains  constant  for  Re  >  20,000.  Additionally,  the 
measurements  revealed  a  “mixing  virtual  origin,”  defined  as  the  far-field  flame 
length  extrapolated  to  </>  =  0,  that  increases  with  increasing  Re  for  Re  <  20,000 
and  then  decreases  with  increasing  Re  for  Re  >  20, 000 .  A  separate  set  of  exper¬ 
iments  indicated  that  the  runs  described  above  were  momentum  dominated  to  the 
farthest  measuring  station  and  that  the  kinetics  of  the  chemiced  reactions  were  feist 
compared  to  the  characteristic  mixing  time.  The  treinsition  of  the  jet  flow  from  a 
momentum-dominated  to  a  buoyeincy-dominated  regime  weis  identified  in  another 
set  of  experiments. 
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CHAPTER  1 


Introduction 


1.1  Background 

The  momentum-dominated,  free  turbulent  jet  is  a  classical  sheair  flow  that  has 
been  experimentally  investigated  by  many  researchers  {e.g.,  Ruden  1933;  Kuthe 
1935;  Hawthorne  et  al.  1948;  Hottel  1953;  Thring  Sz  Newby  1953;  Ricou  &  Spalding 
1961;  Wilson  Sz  Danckwerts  1964;  Becker  et  al.  1967;  Wygneinski  Fiedler  1969; 
Townsend  1976;  Birch  et  al.  1978;  Becker  ^  Yamazaki  1978;  Lockwood  &  Moneib 
1980;  Chen  &:  Rodi  1980;  Dahm  Sz  Dimotakis  1987,  1990;  Dowling  Dimotakis 
1990;  Miller  Sz  Dimotakis  1991).  Although  it  is  perhaps  one  of  the  simplest  of 
flow  geometries,  understanding  the  associated  entrainment  and  mixing  processes 
has  been  an  ongoing  effort  that  is  still  challenging  both  the  experimentalist  and  the 
theorist. 

In  the  context  of  combustion  systems,  sucli  as  in  furnaces,  gas  turbine  engines, 
and  rocket  engines,  it  is  often  desirable  to  have  the  mixing  and  combustion  pro¬ 
cesses  occur  simultaneously  rather  than  sequentially.  The  turbulent  jet  has  seen 
widespread  use  in  such  systems,  providing  a  practical  and  efficient  means  of  mixing 
fuel  and  oxidizer  in  the  combustion  chamber.  Much  of  the  research  in  jet  mixing  has 
been  motivated  by  the  combustion  community  in  the  hope  of  finding  faster,  more 
efficient  ways  of  mixing  reactants  to  provide  more  useful  energy  per  unit  fuel  and  to 
reduce  certain  undesirable  exhaust  products  through  a  better  controlled  combustion 
process  {e.g.,  NOx  in  high  temperature  flames). 


2 


A  common  practice  in  dealing  with  turbulent  flows  has  been  to  assume  that 
above  a  suitably  high  Reynolds  number,  the  properties  of  the  flow  become  indepen¬ 
dent  of  Reynolds  number.  This  assumption  often  extends  to  mixing  and  has  led 
researchers  to  investigate  turbulent  flows  in  the  laboratory  at  much  lower  Reynolds 
numbers  than  those  of  the  actual  devices  utilizing  the  flow.  The  assumption,  or 
hope,  has  been  that  the  experiment  accurately  captures  the  flow  dynamics  of  the 
practical  device,  even  though  the  length  scales,  velocity,  and  fluid  properties  could 
be  very  different. 

Recent  experimental  evidence  has  motivated  a  closer  look  into  the  assumption 
of  Reynolds  number  independence.  Mungal  et  al.  (1985)  found  that  the  amount  of 
product  formation  in  a  chemically  reacting,  two-dimensioned  shear  layer  decreased  as 
Reynolds  number  increased.  In  their  experiment,  product  formation  was  assumed 
to  be  equivalent  to  molecular  mixing,  and  therefore,  the  amotmt  of  molecularly- 
mixed  fluid  changed  with  Reynolds  number.  For  the  current  jet  experiments,  a 
measurable  quantity  that  could  be  related  to  molecular  mixing  was  desired.  Flame 
length  Lf  has  been  shown  to  be  related  to  the  rate  of  molecular  mixing  in  the 
jet  (e.p.,  Broadwell  1982)  and  can  be  measured  experimentally.  Therefore,  Lf , 
considered  to  be  the  axial  distance  required  to  molecvdarly  mix  every  part  of  jet 
fluid  with  reservoir  fluid  to  at  least  the  stoichiometric  mixture  ratio  <(>,  was  chosen 
as  the  quzmtity  to  measure  in  the  current  experiments. 
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1.2  Current  experiments 

The  goal  of  the  present  research  was  to  design  and  construct  a  High  Pressure 
Combustion  Facility  (HPCF)  and  to  conduct  experiments  in  this  facility  to  measure 
average  flame  lengths  in  an  axisymmetric  jet  over  a  wide  range  of  Reynolds  numbers. 
It  should  be  noted  that  there  are  several  methods  for  measuring  average  flame 
lengths,  as  will  be  outlined  in  Ch.  3.  The  partictdar  method  for  determining  L{  in 
these  experiments  was  chosen  not  because  it  is  the  definitive  one  but  because  it  can 
accurately  identify  potentially  small  changes  in  Lf  over  the  large  Reynolds  number 
range  investigated.  That  is,  we  are  primarily  interested  in  determining  if  the  flame 
lengths  change  with  Reynolds  number.  If  there  were  no  Reynolds  number  effects, 
as  conventional  wisdom  would  proclaim,  we  would  expect  a  null  result,  i.e.,  all  the 
measurements  should  collapse  on  a  single  curve  for  the  entire  Reynolds  number 
range  investigated. 

The  jet  coordinate  system  is  shown  in  Fig.  1.1.  The  Reynolds  nximber  of  the 
jet  is  defined  as 

Re  =  £21^  ,  (1.1) 

/^oo 

where  poo  is  the  reservoir  gas  density,  Uo  is  the  jet  exit  mass  flux  velocity  (equal 
to  the  exit  velocity  for  the  current  nozzle),  do  is  the  nozzle  exit  diameter,  and  fi^o 
is  the  reservoir  geis  dyn^lmic  viscosity. 

In  an  effort  to  isolate  Reynolds  number  effects  from  other  effects  (buoyancy, 
kinetics,  etc.),  the  jet  experiment  is  enclosed  in  a  pressure  vessel.  Since  /Zqo  is 
essentially  only  a  function  of  temperature  (e.^.,  Thompson  1984),  the  jet  Reynolds 
number  can  be  varied  by  changing  the  reservoir  gas  density,  which  is  directly  pro¬ 
portional  to  the  reservoir  pressure.  In  this  way,  the  nozzle  exit  diameter  sind  exit 
velocity  can  be  held  fixed,  with  the  jet  in  a  momentum-dominated  regime,  while 
the  Reynolds  number  is  varied  via  pressure. 
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The  facility  is  designed  to  span  2.2  decades  of  pressure,  providing  a  large 
range  of  Reynolds  numbers  for  investigation.  Buoyancy  effects  scale  as 

Ap(x)  ^  Poo-pCj) 

Po  Po 

where  A/?(i)  is  the  average  change  in  density  owing  to  heat  release,  po  is  the  jet 
exit  density,  and  p(x)  is  the  average  density  of  the  jet  cone  up  to  that  x  location. 
These  effects  remain  constzmt  as  the  pressure  is  varied,  because  both  Ap  and  po 
are  linearly  proportional  to  the  pressure.  The  chemiced  kinetic  rate  is  proportional 
to  molecular  number  density  and,  therefore,  for  a  given  concentration  of  reactants, 
the  chemicEil  reaction  rate  gets  faster  as  pressure  (Re)  is  increased. 


5 


Combustion  can  only  occur  when  fuel  and  oxidizer  mix  at  the  molecular  level. 
Therefore,  it  is  of  interest  in  the  present  experiments  to  distingmsh  molecularly- 
mixed  fluid  from  pme  jet  and  reservoir  fluid.  This  is  done  with  the  chemical  reaction 
method  that  was  used  in  water  (acid-base  reaction)  by  Weddell  (in  Hottel  1953), 
Brown  &  Dimotalds  (1976),  Breidenthal  (1981),  Dalim  et  al.  (1984),  and  Koochesfa- 
hani  &  Dimotalcis  (1986),  cind  in  gases  by  MimgaJ  Sc  Dimotahis  (1984),  Hermanson 
Sc  Dimotalcis  (1989),  and  Frieler  Sc  Dimotakis  (1988).  See  Dimotakis  (1989)  for 
a  review  of  these  experiments.  The  method  relies  on  molecular  mixing  to  form  a 
chemical  product  that  can  be  detected  on  a  macroscopic  scale,  e.g.,  heat  release  or 
temperature  rise.  For  the  present  experiments,  the  reactions  are  exothermic,  which 
means  that  chemical  product  formation  is  accompanied  by  heat  release.  If  the  adi¬ 
abatic  flame  temperature  rise  ATf  is  kept  low  enough,  then  the  heat  release  can 
be  expected  to  have  little  effect  on  the  flow  dynamics.  The  adiabatic  flame  tem¬ 
perature  rise  ATf  is  defined  as  the  temperature  rise  achieved  if  all  reactants  burn 
to  completion,  adiabatically,  at  constant  pressure.  For  a  fast  chemical  reaction, 
t.e.,  one  in  which  the  time  lag  between  molecular  mixing  and  product  formation 
is  negligible,  there  is  a  one-to-one  correspondence  between  product  formation  and 
moleculzn  mixing.  In  this  way,  the  heat  release  is  a  tracer  for  molecularly-mixed 
fluid. 

The  advantages  of  the  F2/NO  chemical  system  (described  in  Mungal  Sc  Di¬ 
motakis  1984)  made  it  a  logical  choice  for  the  current  experiments.  This  chemical 
system  is  completely  hypcrgolic.  z.c.,  the  reactants  burn  on  contact  since  they  can¬ 
not  coexist  at  the  molecular  level.  This  obviates  the  need  for  an  ignition  source, 
flameholders,  or  other  flow  disturljing  devices.  Unlike  hydrocarbon  reactions  in¬ 
volving  hundreds  of  complicated  reactions,  the  simple  F2/NO  chemistry  involves 
only  three  intermediate  reactions,  allowing  the  chemical  kinetics  and  the  fluid  me¬ 
chanics  to  remain  tractable.  The  high  kinetic  rate  of  this  system  ensures  that  the 
rate  of  product  formation  is  hmited  by  the  turbulent  mixing  process  and  not  by  the 


chemical  kinetics.  This  system  also  has  the  capability  of  producing  very  low  adia¬ 
batic  flame  temperatures  while  maintaining  fast  chemistry.  This  is  desirable  in  the 
design  of  the  diagnostics,  for  the  reduction  of  buoyancy  effects,  and  in  eliminating 
the  influence  of  finite  kinetic  rate  effects. 

The  F2/NO  chemistry  that  was  employed  in  this  experiment  can  be  described 
by  the  “effective”  reaction: 

2N0 -h  F2  ^  2N0F,  AQ  =  -75-^^^.  (1.3) 

moleF2 

The  reaction  actually  consists  of  two  chain  reactions  with  different  rates  and  heats 
of  reaction,  as  described  in  Sec.  4.3.  The  reactants  provide  a  wide  range  of  heat 
release  while  retzuning  an  overall  reaction  rate  in  the  turbulent  jet  that  is  kinetically 
fast  even  at  the  lower  concentrations. 


1.3  Outline 

Constant  current  resistance  thermometers  were  used  to  mezisure  the  average 
heat  release  in  a  steady,  axisymmetric,  chemically-reacting,  gas-phase,  momentum- 
dominated  jet  discharging  into  a  quiescent  reservoir.  The  measurements  span  the 
downstream  range  from  30  to  240  jet  exit  diameters,  at  Reynolds  numbers  from 
10,000  to  150,000.  Experiments  at  several  stoicliiometric  mixture  ratios  <p  were 
performed  at  each  Reynolds  number  to  identify  possible  Reynolds  number  effects 
on  flame  length.  Chapter  2  describes  the  experimental  facility  and  Ch.  3  describes 
the  measurement  technique.  The  experimental  results  are  discussed  in  Ch.  4,  which 
covers  the  effects  of  Reynolds  number  on  flame  length  zmd  addresses  the  issues  of 
buoyeuicy  and  chemiczd  kinetic  rate. 
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CHAPTER  2 

Experimental  Facility 


2.1  The  main  apparatus 

These  experiments  were  conducted  in  the  newly  constructed  High  Pressure 
Combustion  Facility  (HPCF),  shown  schematically  in  Fig.  2.1  and  pictured  in 
Fig.  2.2.  This  blowdown  facility  was  designed  to  study  turbulent,  reacting  jets 
at  both  reduced  and  elevated  pressures.  In  addition  to  the  capability  of  using  con¬ 
ventional  gaseous  fuels,  the  facility  was  designed  to  handle  more  exotic  reactants 
like  fluorine,  nitric  oxide,  and  hydrogen.  Designing  a  facility  to  safely  use  these 
exotic  gases  presented  many  unique  problems.  Fluorine  is  a  highly  reactive  and 
toxic  gas  that  requires  special  handling  techniques.  Much  of  the  gas  delivery  design 
was  based  on  the  work  of  Mungal  (1983).  Since  the  current  experiments  involved 
higher  concentrations  and  pressures  of  these  gases  than  used  by  Mungal,  fluorine 
compatability  tests  at  high  pressure  euid  high  concentration  were  performed.  These 
materials  tests  were  conducted  with  the  help  of  TRW  at  their  facility  in  Redondo 
Beach. 

The  main  apparatus  is  the  High  Pressure  Reactant  Vessel  (HPRV),  designed  by 
Godfrey  Mimgal  in  collaboration  with  Paul  Dimotakis.  The  HPRV  is  a  107  cm  ( 42 
in.)  internal  diameter  by  183 cm  (72  in.)  length  pressure  vessel  with  an  interior 
volume  of  1.43  m^  (50.5ft^ ).  The  HPRV,  shown  in  Fig.  2.3,  is  capable  of  operating 
at  pressures  from  0.1  atm  to  15  atm  and  is  intemeilly  coated  with  pinhole-free  PFA 
teflon  for  corrosion  resistance.  The  vessel  consists  of  a  shell  with  removable  heads 
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Fig.  2.2  Photograph  of  the  High  Pressure  Combustion  Facility.  The  HPRV  is  in 
the  upper  left  foreground,  with  the  fluorine  reactzmt  tank  directly  beneath 
it.  The  control  panel  used  for  loading  the  jet  eind  reservoir  geises  is  in  the 
right  backgroimd,  and  the  computer  terminal  used  for  data  acquisition 
and  run  sequencing  is  in  the  right  foreground. 


Fig.  2.3  High  Pressure  Reactant  Vessel  schematic.  Height  =  lS3cm  (6  ft),  internal 
diameter  =  107  cm  (3.5  ft). 
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at  each  end  for  full  interior  access.  The  heads  contain  four  7.62  cm  (3  in.)  flanges 
for  possible  coflow  and  one  15.24  cm  (6  in.)  flange  for  jet  nozzles  or  exhaust /opticzJ 
access.  The  shell  contains  six  7.62cm  (3  in.),  diametrically  opposed  flanges  for  pos- 
s.Lie  line-of-sight  absorption  or  ignition  capabilities,  four  3.81cm  (1.5  in.)  flanges, 
and  one  40.64cm  (16  in.)  manhole  for  internal  access.  Optical  access  is  ”ia  three 
25.4cm  (10  in.)  diameter  by  6.35  cm  (2.5  in.)  thick  pyrex  windows  held  in  specially 
designed,  teflon-coated  window  units  manufactured  by  Pressure  Products.  The  tank 
specifications  are  given  in  Appendix  A. 


2.2  The  gas  delivery  system 

The  reservoir  (HPRV)  is  typicedly  loaded  to  the  desired  concentration  of  NO, 
in  a  N2  diluent,  via  intermediate  mixing  vessels  of  known  volume,  utilizing  partial 
pressure  techniques.  The  jet  gas  is  contained  in  a  fluorine  reactant  tank  with  an 
internal  teflon  expulsion  diaphragm  separating  the  jet  charge  from  the  inert  driver 
geis.  A  given  concentration  of  F2  in  a  N2  diluent  is  loaded  into  the  reactEint  tank 
using  the  partial  pressure  technique.  Accuracy  with  the  NO  and  F2  loading  is 
maintained  over  the  large  range  of  pressures  and  concentrations  by  staging  the 
mixing  vessels,  i.e.,  a  small  volume  vessel  is  used  for  sub- atmospheric  conditions  zuid 
a  larger  volume  vessel  is  used  for  atmospheric  and  high  pressure  runs.  The  driver 
side  of  the  expulsion  diaphragm  is  coimected  to  a  nitrogen  surge  tank,  comprised 
of  four  size  lA  cylinders  for  a  combined  internal  volume  of  0.176  m^  (6.2ft^)  of 
driver  g2is.  The  surge  tank  is  pressurized  (typically  35 atm)  to  provide  a  nearly 
constant  stagnation  pressure  to  three  parallel  sonic  metering  valves  of  various  flow 
rates.  The  metered  sonic  N2  stream  supplies  a  fixed  mass  flow  rate  to  the  driver 
side  of  the  expulsion  diaphragm,  resulting  in  a  constant  mass  flow  rate  of  jet  gas 
through  the  nozzle  in  the  HPRV. 
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The  jet  issues  from  the  bottom  of  the  HPRV  upweurds  through  a  2.54  mm  (0.1 
in.)  exit  diameter,  anodized  aluminum  nozzle.  The  nozzle  contraction  section  is 
described  by  a  order  polynomial,  optimized  to  avoid  Taylor- Gortler  vortices 
in  the  concave  parts  emd  for  minimum  exit  boundarj'  layer  thickness.  The  nozzle 
design  method  was  developed  by  Dimotakis  and  is  described  in  Dowling  (1988).  A 
turbulence  management  section  is  located  upstream  of  the  contraction  section.  This 
consists  of  a  3.81  cm  (1.5  in.)  diameter  by  5.08  cm  (2  in.)  length  of  0.318  cm  (0.125 
in.)  core  stainless  steel  honeycomb,  followed  by  four  stainless  steel  screens  of  mesh 
sizes  from  15  to  100  wires  per  inch.  This  results  in  an  exit  turbulence  level  measured 
to  be  less  than  0.3%  over  the  velocity  range  investigated.  The  dyneimic  pressure  is 
set  prior  to  a  run  and  monitored  during  the  run  with  a  Barocel  differential  pressure 
transducer  connected  across  the  nozzle  plenum  and  the  nozzle  exit.  The  nozzle 
is  mounted  on  a  stainless  steel  base  plate  with  three  adjusting  screws  to  provide 
angular  alignment.  The  nozzle  geometric  axis  was  aligned  once  with  a  laser  to  within 
a  few  tenths  of  a  degree  of  the  geometric  centerline  of  the  HPRV  and  secured  to 
prevent  any  subsequent  movement.  Upon  periodic  inspections,  the  nozzle  alignment 
was  found  to  remain  unchanged. 

The  momentum  dizuneter  of  the  jet  d* ,  rather  than  the  nozzle  exit  diameter 
do ,  is  used  to  normalize  the  downstream  coordinate  i .  It  is  defined  by 


d*  = 


^<r/?oo«fo 


where  rho  is  the  nozzle  exit  mass  flux,  Poo  is  the  density  of  the  reservoir  fluid,  and 
Jo  is  the  nozzle  exit  momentum  flux.  This  definition  was  introduced  in  a  limited 
way  by  Thring  &:  Newby  (1953),  used  by  Avery  ic  Faeth  (1974),  and  modified  to  its 
present  form  by  Dahm  &  Dimotakis  (1987).  The  momentum  diameter  can  be  used  to 
collapse  the  results  of  many  different  jet  experiments  (sec  Dahm  i:  Dimotakis  1987) 
aind  is  considered  to  be  the  proper  length  scale  for  normalizing  the  downstream 
distance  in  the  jet.  The  factor  of  2  and  the  y/??  in  Eq.  2.1  were  chosen  so  that  d* 
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reduces  to  the  nozzle  exit  diameter  do  for  density-matched  jet  and  reservoir  fluids 
and  a  perfect  “top- hat”  exit  velocity  profile.  Using  the  jet  and  reservoir  mixture 
densities  and  assiiming  a  uniform  exit  density  profile,  an  axisymmetric  Thwaites 
calculation  was  used  to  estimate  the  boundary  layer  corrections  to  the  velocity 
profile  at  each  rim  condition.  For  the  experiments  of  Sec.  4.1,  the  calculations  give 
d*  =  0.960  do  for  the  Re  =  10,000  runs  and  d*  =  0.991  do  for  the  Re  =  150,000 
runs. 


It  is  also  common  practice  to  allow  for  a  virtual  origin  xo  in  defining  the 
effective  downstream  coordinate  x.  The  virtual  origin  is  determined  by  the  far- 
field  behavior  of  the  jet  and  is  a  correction  for  the  influence  of  the  neeir  field  of 
the  jet.  The  external  shape  of  the  nozzle,  the  jet  exit  turbulence  level,  and  the 
exit  velocity  profile  are  a  few  parameters  that  affect  xo .  Dowhng  (1988)  found  a 
virtual  origin  of  xq  =  0.5  do  for  his  nozzle  at  Re  =  16,000.  The  virtual  origin  for 
the  current  nozzle  was  not  measured.  However,  since  the  nozzle  was  designed  with 
the  same  method  as  Dowling’s,  eind  operates  at  Re  >  10,000,  xo  is  assumed  to 
be  comparable,  t.e.,  Xq  ~  O(do).  It  is  therefore  considered  negligible  zuid  will  be 
ignored  throughout  this  study. 


2.3  The  pressure  relief  and  exhaust  systems 

The  facility  operates  in  a  blowdown  mode  under  constant  pressure  conditions. 
Pressure  relief  for  the  HPRV  is  provided  by  a  3.81  cm  (1.5  in.)  Camflex  relief  valve 
equipped  with  a  Tajdor  Microscan  500  Controller.  The  Controller  uses  a  feedback 
control  system  that  monitors  the  pressure  in  the  HPRV,  compares  it  to  a  set-point 
pressure,  auid  opens  the  Camflex  in  a  preprogrammed  manner  to  maintain  the  set- 
point  pressure.  The  Camflex  is  used  only  for  run  pressures  >  latm,  while  for 
sub- atmospheric  runs,  the  Camflex  remains  shut  and  an  acceptable  pressure  rise 
occurs  in  the  HPRV.  The  exhaust  from  the  Camflex  is  plumbed  into  a  downstream 


catch  bag  that  is  shared  with  the  HF  Shear  Layer  Combustion  Facility.  The  catch 
bag  (described  in  Hall  1991)  is  equipped  with  NaOH  and  water  showers  to  neutralize 
the  exhaust  products. 

No  commercially  available  pressure  transducer  had  the  dynzimic  range  to  sense 
the  pressure  rise  in  the  HPRV  over  the  entire  relief  range  (latm  to  15 atm)  with 
the  resolution  required.  It  was  therefore  necessary  to  create  a  “sliding”  pressure 
measurement  that  had  verj'  high  resolution  over  a  smeill  dynamic  range,  but  whose 
range  could  be  centered  anywhere  within  the  operating  pressmre  of  the  HPRV.  This 
was  achieved  by  storing  the  HPRV  starting  pressure  in  a  thermally-insulated  volume 
that  was  large  compared  to  the  internal  volume  of  a  Barocel  differential  pressure 
transducer.  This  volume  was  ported  to  one  side  of  the  differential  transducer  and 
isolated  via  a  valve  before  the  run  began.  The  other  side  of  the  transducer  was 
connected  to  the  HPRV  resulting  in  a  highly  accurate  (0.006  psi  resolution)  pressure 
measurement  that  could  “slide”  anywhere  within  the  pressure  range  of  the  HPRV. 
This  arrangement  performed  well,  with  the  Camflex  providing  constant  pressure  in 
the  HPRV  to  within  0.5%,  or  better,  of  the  rim  pressures  investigated  ( 1  atm  to 
15  atm). 

2.4  Diagnostics 

« 

The  diagnostics  utilized  in  these  experiments  were  designed  to  give  an  accurate 
measure  of  average  flame  length.  Temperature  measurements  were  performed  with 
a  set  of  16  long  cold  wires  opiated  as  constant  current  resistance  thermometers. 
Each  wire  consists  of  a  25 /tm  (0.001  in.)  diameter  by  57.8cm  (22.75  in.)  length 
of  Pt-10%Rh  wire  welded  to  a  pair  of  356 //m  (0.014  in.)  diameter  by  5cm  (2  in.) 
long  inconel  prongs.  Platinum  was  chosen  for  its  high  coefficient  of  resistivity  and 
rhodium  was  chosen  to  increase  the  strength.  The  prongs  serve  as  cantilever  springs 
to  maintain  tension  in  the  wires  and  also  to  provide  the  electric^J  connection  to  the 
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Fig.  2.4  Schematic  of  the  diagnostics  with  the  coordinate  axes  indicated. 


leads  going  to  the  eunplifier  circuits.  The  prongs  are  attached  to  circular  teflon 
collars  mounted  on  two  vertical  stainless  steel  rods  located  inside  the  HPRV.  The 
wires  are  stretched  across  the  jet  centerline  from  a-/rfo  =  30  to  240  and  are  spaced  at 
axial  locations  selected  to  give  equal  spacing  on  a  logio(x/do)  axis.  This  spacing  was 
suggested  by  Paul  Dimotakis  in  view  of  his  conjecture  of  a  logarithmic  behavior  of 
the  temperature  mezisurements  in  the  fleime  region  (private  communication).  Figure 
2.4  shows  a  schematic  of  the  wires  with  the  coordinate  axes  indicated. 


Each  wire  has  a  typical  resistance  of  230  ohms  and  is  provided  with  a  constant 
current  of  0.2  mA  by  an  zunplifier  circuit  originally  designed  by  Dan  Lang  and  Paul 
Dimotalds,  zind  modified  by  Cliff  Frieler.  Under  natural  convection,  a  worst  case, 
this  yields  an  ohmic  resistance  heating  temperature  rise  of  0.0006  K.  Therefore, 
the  ohmic  heating  is  negligible  when  compared  to  the  temperatures  measured  in 
the  experiments  zmd  the  wire  sensitivity  is  velocity  independent.  Because  of  the 
extremely  large  aspect  ratio  of  the  wires  (22,750)  zmd  the  fact  that  the  prongs  are 
located  weU  outside  of  the  jet  cone,  end  effects  are  also  negligible. 

The  circuit  measures  wire  resistance  changes  due  to  heating  by  the  reacting 
jet  flow,  and  then  amplifies  the  signal  via  a  variable  gziin  stage.  The  signzil  is  then 
filtered  with  a  third-order  low-pass  Butterworth  filter  with  the  knee  frequency  set 
at  200  Hz  (see  Table  A. 3  for  estimates  of  the  wire  frequency  response).  The  rms 
noise  without  flow  was  measured  at  ~  0.015  K  and,  with  the  long  averaging  times 
involved,  did  not  affect  the  mean  temperature  measurements.  The  measurement 
technique  is  further  described  in  Ch.  3. 


2.5  Data  acquisition 

Data  acquisition  was  handled  by  a  DEC  LSI-11/73  CPU-based,  RT-11  com¬ 
puter  system  together  with  a  Data  Translation  DT3362  12-bit,  32  differential  input 
channel  A/D  converter  with  a  maximum  data  rate  of  250  kHz.  The  system  also 
has  an  Andromeda  Systems  DAC-11  D/A  converter,  used  for  operating  the  remote 
valves,  and  a  Wavetek  Model  27S  fimction  generator,  whicli  provided  the  clock 
signal  for  the  A/D  board. 

The  16  cold  wires,  two  reference  cold  wires,  the  nozzle  dynamic  pressure,  and 
the  HPRV  pressure  were  sampled  at  500  Hz  each  for  a  combined  data  rate  of  10  kHz. 
A  typical  run  time  of  20  seconds  resulted  in  200,000  individuzd  measurements.  The 
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computer  also  performed  the  valve  sequencing  during  a  run.  This  consisted  of 
starting  the  A/D  data  acquisition,  opening  a  firing  valve  to  initiate  the  sonic  N2 
stream,  waiting  a  preset  amount  of  time,  and  then  opening  the  nozzle  valve  to  start 
the  experiment.  After  a  designated  run  time  elapsed,  the  computer  simultaneously 
closed  all  valves  and  stopped  the  A/D  data  acquisition.  After  a  run,  the  data  were 
transferred  via  Ethernet  to  the  HYDRA  VAX/ VMS  cluster  for  processing  (described 
in  Appendix  B). 


CHAPTER  3 


Measurement  technique  and  jet  scaling 


Various  methods  have  been  used  to  measure  flame  lengths.  For  exeimple,  Dahm 
ti  al.  (1984)  used  planar  laser-induced  fluorescence  with  an  acid-base  reaction  in 
water.  They  seeded  the  basic  jet  fluid  with  a  dye  that  would  fluoresce  when  excited 
by  a  lEiser  sheet  imtil  the  jet  dye  was  mixed  to  a  chosen  stoichiometric  mixture  ratio 
with  acidic  reservoir  fluid,  whereupon  the  fluorescence  would  “turn  off”.  By  slicing 
the  jet  centerline  with  a  laser  sheet  oriented  parallel  to  the  jet  Jixis,  they  were  able 
to  visuzdize  a  two-dimensional  plane  containing  the  axis  of  the  jet.  They  meaisured 
the  farthest  downstream  location  where  fluorescing  dye  could  be  seen  in  the  laser 
sheet  and  called  this  the  instantaneous  flame  length.  They  then  ensemble- averaged 
many  of  these  instantaneous  measurements  to  define  an  average  flame  length. 

A  second  method  is  to  locate  the  farthest  downstream  location  where  the  ra¬ 
diated  light  from  high  temperature  soot  formed  in  hydrocarbon  flames  cEin  be  seen 
{e.g.,  Zukoski  et  al.  1984;  Toner  1987).  Soot  is  formed  bls  a  by-product  of  the  com¬ 
bustion  process,  and  radiates  light  only  when  sufficiently  hot.  It  is  aissumed  that 
the  soot  radiates  light  only  within  the  flame  zone,  where  the  highest  temperatures 
exist,  and  that  beyond  the  flame  tip,  the  soot  cools  to  a  temperatiure  where  visible 
radiation  ceases.  Therefore,  the  last  occunence  of  visible  soot  should  mark  the  end 
of  the  flame  tip.  There  is  some  concern  over  whether  the  soot  cools  immediately 
beyond  the  flame  tip,  or  whether  it  continues  to  radiate  light  for  some  distance 
beyond  the  flame  tip  as  it  is  cooling.  The  latter  condition  could  result  in  over¬ 
estimated  flame  lengths.  Additional!}',  high  heat  release  flzunes  are  also  typically 
buoyancy-dominated . 
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Others  have  judged  visued  average  flame  lengths  from  the  radiated  light  of 
a  flame  by  eye  or  used  long  exposure  photographs  of  the  flame  against  a  black 
background  {e.g.,  Hawthorne  ei  al.  1948;  Becker  k.  Liang  1978). 

In  the  current  experiments,  we  accurately  measured  the  low  heat  release  of  the 
chemical  reaction  to  estimate  flame  lengths.  The  diagnostics  utilized  in  the  present 
experiments  are  designed  to  give  a  sensible,  accurate,  and  consistent  measure  of 
flame  length  using  constant  current  resistance  thermometers. 

As  noted  in  Sec.  2.4,  the  thermometers  are  long,  thin  Pt-10%Rh  cold  wires 
stretched  across  the  jet  centerline  at  16  axial  locations.  If  ohmic  heating  of  the 
wires  and  the  slight  nonlinearities  in  the  temperature  coefficient  of  resistance  are 
neglected,  then  the  wire  voltage  chzinge  AV  for  a  given  temperature  change  AT 
can  be  given  as 

AV  =  IRoaAT  ,  (3.1) 


where 

I  =  current  through  the  wire, 
a  =  temperature  coefficient  of  resisteince, 

Ro  =  probe  cold  resistance. 

The  wire  resistance  R  is  well  represented  as  a  linear  fimction  of  temperature,  given 
by  i?  =  i2o(l  +  qAT).  See  Sandborn  (1972).  This  zissumption  is  valid  for  the 
temperature  range  of  the  current  experiments,  and  the  voltage  increase  is  therefore 
lineeu-  with  temperature. 


X 


AT(x,y) 


Fig.  3.1  Schematic  of  the  jet  flame  zone. 


3.1  Jet  Scaling 

The  rezison  for  choosing  this  particular  measurement  technique  was  motivated 
by  the  scaling  laws  for  an  axisymmetric  free  jet.  The  resistance  of  the  cold  wires, 
when  stretched  across  the  flame,  yields  an  accurate  estimate  of  the  line-integrated 
heat  release,  i.e.,  temperature  rise,  resulting  from  the  chemical  reaction.  At  the 
end  of  the  flame  zone,  shown  schematically  in  Fig.  3.1,  the  fuel  hzis  been  completely 
consumed,  so  no  heat  release  occurs  downstream  of  this  location.  Beyond  the  flzune 
tip,  the  existing  heat  is  simply  diluted  by  the  subsequently  entrained  reservoir  fluid. 
Therefore,  for  x  >  Li  (where  x  is  the  streamwise  (jet  axis)  coordinate  and  Lf  is 
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the  average  flame  length),  the  temperature  excess  is  a  conserved  scalar.  Assuming 
that  the  flow  achieves  a  self-similar  temperature  profile  beyond  the  flame  tip  (as 
observed  by  Becker  Yamazalci  1978),  the  temperature  excess  for  x  >  L{  should 
obey  the  jet  similarity  law  for  a  scalar  quantity  in  a  momentum-dominated  jet.  This 
similarity  law  can  be  derived  from  dimensional  arguments  (e.^.,  Chen  &  Rodi  1980), 
has  been  verified  experiment zJly  ( e.^.,  Dowling  ^  Dimotakis  1990),  and  is  given  by 


AT(x,ij) 

ATf 


(3.2a) 


where  AT(i,y)  is  the  temperature  rise  at  (r,  y),  y  is  the  transverse  coordinate, 
ATf  is  the  adiabatic  flame  temperature  rise,  k  is  a  constant  (determined  by  exper¬ 
iment),  d*  is  the  momentum  diameter  of  the  jet  (see  Eq.  2.1),  g{T})  is  the  similarity 
mean  profile  function  (determined  by  experiment),  rj  is  the  similarity  variable  de¬ 
fined  as 

^  -  y 

V  —  ~  » 

X  —  X0 

with  Xo  the  virtual  origin.  The  function  g(7j)  has  a  maximum  value  of  1  at  r;  =  0 
and  goes  to  zero  outside  the  jet  edges  (  t;  >  0.25 ). 


Taking  the  line  integral  of  Eq.  3.2a  along  y,  at  constant  x  (along  a  wire),  and 
realizing  that  y(r/)  goes  to  zero  before  y  =  ±Xw/2,  yields 


1  AT(a-,!,)  ,  d'/U  ,  , 

T  /  — AT — “i'  “  “z — 

■^w  J-Lw/2  X  —  Xo 


d*/L^ 

X  —  Xo 


i: 

-(£-\  r 

\  Ly,  J 


9{v)dy 


g{v)dT] 


7^  fn(.T)  , 


(3.3) 
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where  Lw  is  the  span  (length)  of  the  wire.  The  temperature  measurement  performed 
by  each  wire  is  the  line  integral  at  constant  x  given  by 


U-  J-I 


AT{x,y)dy  . 


This  measurement  is  analogous  to  the  product  thickness  used  for  the  shear  layer 
(e.g.,  Dimotakis  1989).  The  product  thickness  ^p(x),  normalized  by  Xw,  will  be 


defined  as 


Sp(x)  _  J. 

Lv,  L\ 


--f 

iw  J-, 


AT(x,y) 


The  product  thickness  6p(x)  can  be  thought  of  as  the  width  of  a  top-hat  temper¬ 
ature  profile  of  height  ATf  that  results  in  the  same  area  as  that  found  under  the 
AT(x,  y)  profile  for  that  x  location.  Equations  3.5  and  3.3  axe  equivalent  to  Eq.  3.4 
divided  by  ATf.  Therefore,  the  time-averaged,  line-integrated  mean  temperature 
rise  measured  at  each  wire  beyond  the  flame  tip  at  L{  is  expected  to  asymptote  to 
a  constant  \'alue.  Accordingly,  the  line  integral  will  rise  to  that  asymptotic  vadue 
as  dictated  by  the  cumulative  mixing  and  chemical  product  formation  up  to  that 
location  for  values  of  x-  <  L{ . 


Figure  3.2  shows  a  sample  of  the  raw  data  from  the  16  wires  spanning  x/do  = 
30  to  240,  along  with  two  reference  wires  placed  outside  of  the  jet  cone  at  x/do  =  0 
and  240,  amd  the  scaled  velocity  trace.  Starting  with  the  bottom  trace  and  working 
up,  the  lowest  trace  is  the  reference  wire  at  x/do  =  0  (partially  obscured  by  the 
x-aixis),  then  the  16  measurements  from  x/do  =  30  to  240,  the  other  reference 
wire  at  x/do  =  240,  and.  finally,  the  scaled  velocity  trace.  The  reference  wire  at 
i/do  =  0  meaisures  the  temperature  of  the  ambient  reservoir  fluid  and  the  reference 
wire  trace  at  x/do  =  240  indicates  the  arrival  of  hot  recirculated  products  at  the 
Izist  measuring  station,  thus  determining  the  temporal  extent  of  uncontaminated 
data.  The  scaled  velocity  trace  indicates  the  time  required  to  establish  steady  state 
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0  4  8  12  16  20 


t  (sec) 

Each  channel  offset  by  0.75 K 

Fig.  3.2  Sample  traces  of  cold  wire  raw  data  (see  text  for  description  of  plot).  The 
dashed  lines  mark  the  averaging  interval  used  for  Fig.  3.3. 


These  data  were  processed  to  produce  Fig.  3.3.  It  plots  the  nondimensionaJ 
product  thickness,  given  by  Eq.  3.5,  versus  logjo(x/d*)  for  four  different  <^’s  at 
ATf  =  7K.  The  solid  line  for  each  set  of  data  is  a  compound  curve  fit  with  a 
linear  leeist  squares  fit  to  the  ramp  region  and  a  cubic  spline  fit  to  the  knee  region, 
with  the  appropriate  matcliing  of  slope  and  curvature  at  the  crossover.  The  fieime 
length  L{  will  be  defined  eis  the  99%  of  maximum  6p/£w  point,  in  the  same  spirit 
as  the  definition  of  a  99%  boundary  layer  thickness.  The  expected  increase  in  the 
integrated  temperature  rise  is  observed,  as  well  as  the  aisymptotic  value  at  the  end 
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logjoCx/d*) 


Fig.  3.3  Sample  plot  of  product  thickness  versus  logio(x/d*)  for  severeil  (t>'s. 


of  the  fleime  tip.  The  measurements  display  a  linear  behavior  in  the  ramp  region 
when  plotted  in  the  semi-logarithmic  coordinates  suggested  by  the  choice  of  wire 
spacing,  and  corroborate  the  conjecture  by  Dimotakis  mentioned  in  Sec.  2.4. 


The  fact  that  the  meeisurements  achieve  a  constant  level  as  predicted  suggests 
that  the  sissumption  of  a  self-similar  temperature  profile  beyond  the  flame  tip  is  ® 

v2ilid.  It  should  be  noted  that  based  on  the  current  definition  of  L( ,  1%  of  the 
heat  release  occurs  beyond  x  =  L[ .  It  is  possible  that  this  additiouEil  heat  release 
could  conspire  with  a  non-  .milar  temperature  profile  that  is  developing  into  the  self-  • 

similar  one  to  result  in  a  constemt  temper:  ture  measurement.  This  might  cause  us  to 
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underestimate,  or  overestimate,  the  flame  length,  as  defined  by  some  other  criterion. 

,  However,  we  axe  mainly  interested  in  changes  in  Li  with  Reynolds  number,  and 

therefore,  regardless  of  the  deflnition  chosen,  it  is  the  change  in  Li  with  increasing 
Reynolds  number  that  is  important  here. 

*  It  should  also  be  noted  that  the  constant  level  achieved  by  the  temperature 

measurements  beyond  the  flame  tip  in  Fig.  3.3  only  apphes  to  a  momentum- 
dominated  jet.  The  entrainment  rate  of  a  buoyant  jet  is  higher  than  that  of  a 
►  momentmn-dominated  jet  {e.g..  Becker  Yamazaki  1978).  This  additional  en¬ 

trainment  of  cold  reservoir  fluid  will  typically  cause  the  temperature  measurements 
in  the  buoyancy-dominated  region  to  decreeise  rather  than  remain  constant.  The 
^  fact  that  the  measureme.its  in  Fig.  3.3  achieve  a  constant  level  demonstrates  that 

this  choice  of  flame  temperature  places  the  jet  in  the  momentum-dominated  regime 
out  to  the  farthest  measuring  station  for  all  the  (f)'s  investigated.  This  will  be 
^  discussed  further  in  Sec.  4.2. 


► 

) 

) 
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CHAPTER  4 

Results  and  Discussion 


The  present  experiments  investigate  a  steady,  effectively  unconfined,  turbulent 
jet  issuing  into  a  quiescent  reservoir.  One  jet  nozzle  with  an  exit  diameter  of  do  = 
2.5  mm  (0.1  in.)  was  used  for  all  of  the  experiments  reported  here.  For  reasons 
discussed  in  Sec.  4.3,  all  runs  were  performed  at  pressures  between  latm  <  p  < 
15 atm.  The  jet  g£is  was  always  F2  in  a  N2  diluent  and  the  reservoir  gas  was  sJways 
NO  in  a  N2  diluent. 


4.1  Reynolds  Number  Effects 


Experiments  at  four  stoichiometric  mixture  ratios  4>  were  performed  at  each 
Reynolds  number  to  produce  four  different  flame  lengths.  The  stEmdeud  diffusion 
flame  convention  {e.g.,  Kuo  1986)  defines  <f>  as  the  fuel  to  oxidizer  ratio,  divided  by 
the  fuel  to  oxidizer  stoichiometric  ratio,  i.e.. 


6  = 


(4.1a) 


where  /  Euid  o  axe  mass  or  mole  fractions.  There  is  no  need  to  distinguish  which 
species  is  fuel  and  which  is  oxidizer  for  this  study.  But,  in  an  effort  to  be  consistent 
with  diffusion  flame  theory  in  which  the  fuel  is  typically  carried  by  the  jet  gas 
and  the  quiescent  reservoir  contains  the  oxidizer,  tliis  work  will  follow  the  same 
convention.  Therefore,  <i>  will  be  defined  as  the  mass  fraction  of  F2  in  N2 ,  divided 
by  the  mass  fraction  of  NO  in  N2 ,  normalized  by  the  corresponding  stoichiometric 
mass  fraction  ratio,  or 

.  ^  (Ff^/Ino) 

"  (rF,/rNo)st  ’ 


(4.1b) 
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logjo(x/d*) 

Fig.  4.1  Normalized  product  thickness  versus  logio(ar/d*),  <^  =  7. 


where  the  meiss  fraction,  i'i ,  is  defined  for  a  mixture  of  N  species  as 


= 


m; 


N 

£ 

1=1 


(4.2) 


2md 


rrii  =  mass  of  the  i’’*' species  . 


The  adiabatic  flame  temperature  rise  ATf  for  each  <j>  was  computed  using  the 
CHEMKIN  chemical  kinetics  progrmn  (Kee  et  al.  1980).  The  four  values  of  <i>  with 
the  corresponding  meiss  fractions  and  ATf ’s  are  given  in  Table  4.1. 
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Fig.  4.2  Normzilized  product  thickness  versus  logjo(x/d*),  4> 


Table  4.1:  Stoichiometric  mixture  ratios,  mass  fractions,  and 
adiabatic  flame  temperature  rises  used  in  Sec.  4.1. 


NO  ATf  (K) 


0.0029 


0.0028 


0.0027 


0.0027 


<i> 

7.0 

0.0127 

10.0 

0.0174 

14.0 

0.0237 

17.9 

0.0300 

The  choice  of  A  T{  wjis  determined  from  several  preliminary  experiments.  The 
lowest  value  of  ATf  that  could  be  measured  accurately  and  with  repeatability  was 
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Fig.  4.3  Normsdized  product  thickness  versus  logi0(i/d*),  ^  =  14  and  15. 


sought  to  reduce  both  buoyancy  effects  and  heat  release  effects  on  the  flow  dynamics. 
Section  4.2  describes  the  choice  of  ATf  «  7K  on  the  basis  of  buoyancy  concerns 
and  Sec.  4.3  demonstrates  that  this  low  flame  temperature  is  characterized  by  a  feist 
kinetic  rate. 

Figure  4.1  shows  the  data  for  the  4>  =  7  nms.  The  solid  line  is  a  compound 
curve  flt  with  a  lineeir  least  squares  fit  in  the  ramp  region  and  a  cubic  spline  fit  to 
the  knee  region.  The  run  at  Re  =  10, 000  required  a  higher  adiabatic  flame  tem¬ 
perature  because  of  fluorine  depletion  in  the  nozzle  (described  further  in  Sec.  4.3). 
Similar  results,  for  the  runs  at  0  =  10,  14,  and  18,  are  given  in  Figs.  4.2,  4.3,  and 


Fig.  4.4  Normalized  product  thickness  versus  logio(x/d*),  <t>  =  18. 


4.4,  respectively.  The  runs  were  fairly  repeatable  as  shown  by  the  good  agreement 
of  the  repeated  runs  at  Rt  =  20,000  (squares)  in  Fig.  4.2.  All  of  the  runs  exhibit 
a  linear  ramp  region  when  plotted  in  semi-logarithmic  coordinates  eind  a  neeirly 
const£mt  zisymptotic  level  after  the  knee  or  flame  tip  region.  Figure  4.5  shows  the 
entire  set  of  data  from  Figs.  4.1,  4.2,  4.3,  and  4.4,  plotted  on  the  same  scale. 


The  average  fleime  lengths  were  estimated  from  these  curves  by  determining 
the  99%  of  maximum  point.  These  flame  lengths  are  plotted  versus  <f>  for 

the  various  Reynolds  numbers  in  Fig.  4.6.  The  straight  lines  are  linear  least-squares 
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Fig.  4.5  Normalized  product  thickness  versus  logio(x/d*)  for  the  various  <^’s. 


fits  to  the  data  of  the  form 

Lt 

—  —  A(i>-\-  Lq  ,  (4.3) 

where  A  and  Lq  are  functions  of  Re  with  A  being  the  slope  and  Lq  being  defined 
as  the  “mixing  virtual  origin.”  This  mixing  virtuEil  origin  is  a  new  observation  and 
corresponds  to  the  abscissa-intercept  in  Fig.  4.6,  t.e.,  the  fiame  length  extrapolated 
to  =  0.  If  the  fiow  could  be  described  by  its  fax-field  behavior  at  the  nozzle  exit, 
then  Lq  would  be  the  fiame  length  measured  in  the  limit  oi  <i>  — >  0.  The  variation 
of  Lq  in  I-’ig.  4.6  indicates  that  this  mixing  virtual  origin  systematically  decreases 
with  increasing  Reynolds  number  for  Re  >  20,000.  Some  of  the  variation  in  Lq 


Fig.  4.6  Normalized  flame  length  versus  <i>. 


could  be  attributable  to  the  uncertainty  in  the  curve  fits  used  for  the  flame  length 
determination.  Another  possibility  could  be  a  variation  in  iq  with  increasing  Re . 
However,  even  combined,  these  effects  would  not  be  leirge  enough  to  explain  such  a 
large  variation  in  I-o  •  A  nonzero  intercept  in  itself  is  not  surprising,  if  one  rezdizes 
that  the  potential  core  extends  roughly  six  diameters  from  the  nozzle  exit  and  the 
self-similzir  region  of  the  jet  typically  begins  around  x/do  =  20  (e.^.,  Dowling 
&:  Dimotakis  1990).  Therefore,  fully  developed  turbulent  mixing  on  a  moleculzu- 
scale  requires  some  distance  downstream  of  the  nozzle  exit  to  develop,  hence  the 
streamwise  delay  in  moleculzir  mixing.  The  fact  that  Lo  hzis  such  a  large  value  and 
that  it  varies  with  Re ,  however,  is  a  new  observation. 


Lf/d 
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Fig.  4.7  Comparison  of  normalized  fleime  lengths  versus  <j>. 


The  data  of  Fig.  4.6  are  compared  with  the  average  flame  length  measurements 
of  Weddell  (in  Hottel  1953),  Dahm  et  al.  (1984),  and  Dowling  (1988)  in  Fig.  4.7. 
The  measurements  of  Weddell  and  Dahm  et  al.  were  performed  in  water  with  acid- 
base  reactions,  utilizing  method.s  where  the  jet  fluid  could  be  seen  until  mixed  to 
the  stoichiometric  mixture  ratio  with  acidic  reservoir  fluid.  Their  data  indicate 
the  farthest  downstream  location  where  mireacted  jet  fluid  could  be  seen.  DowHng 
worked  in  gas  phase  and  inferred  maximum  fleime  lengths  from  the  probability  den¬ 
sity  functions  (PDFs)  of  his  Rayleigh  scattering  concentration  measurements.  He 
assumed  a  99%  complete  reaction  in  his  ceJculations  of  the  flame  lengths  presented 
in  Fig.  4.7. 


The  current  results  yield  mean  flame  lengths  that  are  roughly  25%  shorter 
than  those  of  Dahm  et  al.  and  Dowling  at  the  higher  <j>'s.  This  discrepancy  could 
be  attributed  to  two  possibilities.  The  first  is  the  difference  in  Schmidt  number, 
since,  for  water.  Sc  w  10^ ,  and  for  a  gas,  Sc  a;  1 .  However,  the  data  of  Dahm  et 
al.  in  water  agree  with  the  inferred  gas-phase  flame  lengths  of  Dowling,  implying 
that  there  is  no  Sc  effect.  The  other  possibility  could  be  the  different  methods 
used  for  estimating  the  flame  length.  Dalim  et  al.  ensemble-averaged  the  farthest 
downstream  location  where  fluorescing  dye  could  be  seen  to  define  an  average  flame 
length,  and  Dowling  estimated  maximum  flame  lengths. 

An  important  feature  of  Dowling’s  data  is  the  reduction  in  the  slope  A  of 
the  Ll/d*  curve  from  >1  «  9  for  Re  =  5,000  to  A  ss  8.35  for  Re  =  16,000. 
Although  this  change  is  not  zis  large  as  in  the  current  experiments,  within  a  similar 
Reynolds  number  range,  the  behavior  is  qualitatively  simileir. 

The  main  result  of  this  investigation  is  that  there  are  Reynolds  number  effects 
on  flame  length  as  illustrated  in  Figs.  4.S  and  4.9.  The  values  of  the  slope  A  axe 
plotted  in  Fig.  4.8.  A  least-squares  fit  to  the  data  of  Weddell  and  Dahm  et  al.  in 
Fig.  4.7  results  in  a  value  of  -4  rj  9  for  both  sets  of  data.  The  bars  on  the  data 
indicate  the  statistical  confidence  in  the  linear  least-squares  fit.  The  slope  for  the 
data  of  Dowling  changes  roughly  7%  for  a  change  in  Reynolds  number  from  5,000 
to  16,000.  The  data  of  the  current  experiments  display  a  larger  change  in  A,  i.e., 
roughly  13%  from  Re  =  10,000  to  Re  =  20,000.  The  slope  for  the  data  of  the 
current  experiments  reaches  a  constant  value  of  A  rj  5.6  for  Re  >  20,000.  The 
entrainment  studies  of  Ricou  Spalding  (1961)  found  that  the  entrainment  rate  of 
their  gas  jet  became  constant  for  Re  >  25,000.  This  Reynolds  number  is  similar  to 
the  one  where  A  becomes  constant,  as  described  above.  However,  the  behavior  of 
the  mixing  virtual  origin  does  not  seem  to  have  reached  an  asymptotic  behavior  at 
this  Reynolds  number  or  even  up  to  Rt  =  150,000.  Whether  or  not  the  threshold 
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Fig.  4.8  The  slope  A  of  the  curve  fits  to  the  data  of  Fig.  4.7  versus  logjp  Rt.  Beits 
represent  the  statistical  confidence  in  the  linear  leeist-squEires  fit  and  are 
omitted  if  smEiUer  than  the  plotted  symbol. 

Reynolds  number  of  Ricou  ^z  Spalding  can  be  correlated  to  the  asymptotic  behavior 
of  A  is  imclear. 

The  vEuriation  in  Lq  is  plotted  in  Fig.  4.9.  All  of  the  experiments  have  a 
nonzero  mixing  virtual  origin.  The  trend  seems  to  be  for  this  quEuitity  to  increEise 
with  increEising  Reynolds  number  up  to  Re  =  20, 000 ,  Euid  then  to  decrease.  Unfor- 
tvmately,  rims  at  higher  Reynolds  numbers  are  required  to  determine  if  this  quantity 
asymptotes  to  some  constant  le-v'el  with  increasing  Reynolds  number.  It  is  interest¬ 
ing  to  note  that  Lq  accounts  for  an  increEisingly  larger  fraction  of  the  fiEime  length 
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Fig.  4.9  Mixing  virtual  origin  Lo  of  the  curve  fits  to  the  data  of  Fig.  4.7  versus 
logio-Re. 

as  the  flame  length  decreeises,  i.e.,  zis  0  decreases.  Therefore,  changes  in  Lq  with 
Reynolds  number  have  a  proportionally  larger  effect  on  Lf  at  smaller  ^’s. 
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4.2  Buoyancy  Effects 

Buoyancy  is  an  issue  that  weis  explicitlj'^  addressed  in  these  experiments.  The 
similarity  law  of  Eq.  3.2a  is  valid  only  for  the  momentum-dominated  region  of  the 
jet.  The  positively  buoyant  jet  (t.e.,  a  jet  in  which  the  momentum  flux  increases 
with  x)  entrains  at  a  higher  rate  than  the  momentum-dominated  jet  (e.g.,  Becker 
&  Yamazaki  1978).  Therefore,  the  similarity  law  that  results  in  a  const£int  average 
temperature  measurement  beyond  the  flame  tip  will  change.  This  cam  be  imderstood 
by  revisiting  the  scaling  arguments  of  Sec.  3.1. 


Recall  that  the  temperature  excess  beyond  the  flame  tip  is  a  conserved  scalzu 
obeying  the  similarity  law  for  a  momentum-dominated  jet  given  by  Eq.  3.2a.  The 
temperature  excess  beyond  the  flame  tip  for  a  buoyaincy-dominated  jet  is  also  a 
conserved  scalar  and  obeys  the  similarity  law  {e.g.,  Chen  Sz  Rodi  1980)  given  by 

=  »■, - ,  (4-4) 

ATf  (x  — xo)*/'* 

where  0  for  the  buoyant  jet  is  different  from  the  k  of  the  momentum-dominated 
jet  in  Eq.  3.2a  and  the  mean  temperature  profile  fimction  g{rf)  is  also  different. 
The  similarity  variable  ry  is  the  same  as  for  the  momentum-dominated  jet,  given  by 
T)  =  y/{x  —  xq).  Integrating  Eq.  4.4  along  y,  at  constant  x,  and  again  realizing 
that  g{T])  goes  to  zero  before  y  =  ±Lw  /2,  yields 


1  AT{x,y) 


-Lw 


dy  = 


Ly,{x-Xo)^^^  J-L^/2 


/Lyf  /2 

9{v)dy 

■Iw/2 


0d* 

Ly,(x  -  xo)®/^ 


f  g{v)dy 

•/— oo 


Bd*  r°° 

TT^ — ^  /  s{v)dv 

Ly,{x  -xoYi^ 


~  .T 


(4.5) 


logjQ(x/d*) 


Fig.  4.10  Scan  of  ATf  for  <l>  =  17.9,  Re  =  20,000,  Uq  =  62.4m/s,  p  =  2atm. 


Therefore,  the  average  temperature  measurement  will  decay  beyond  the  flame  tip  at 
some  point  when  buoyancy  forces  become  comparable  to  the  jet  source  momentum 
flux. 


For  the  experiments  of  Sec.  4.1,  ATf  was  chosen  by  the  requirement  that  the  jet 
remain  momentum-dominated  up  to  the  Isist  measvu-ing  station  at  i/do  =  240  for 
the  4> ’s  investigated.  This  was  done  by  performing  a  set  of  preliminary  experiments 
at  several  <^’s.  The  longest  flzime  length  occurs  at  0  =  17.9.  This  case  was  exam¬ 
ined  by  systematically  increasing  ATf  until  the  temperature  measurement  beyond 
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logiQ(x/d*) 

Fig.  4.11  Scan  of  ATf  for  4>  =  10,  i?e  =  20,000,  Uq  =  62.4m/s,  p  =  2atm. 


the  flame  tip  was  no  longer  constant,  but  was  observed  to  decrease.  Figure  4.10 
shows  these  4>  =  17.9  runs.  For  the  run  with  ATf  =  10  K,  the  temperature 
measurement  edter  the  flame  tip  never  attains  a  constant  value  but  decreases  with 
increeising  x/d* . 

Figures  4.11  and  4.12  show  similar  flame  temperature  scans  at  <^  =  10  and  7, 
respectively.  These  shorter  flames  are  not  affected  by  buoyancy  until  much  higher 
flame  temperatures.  As  expected,  the  longest  flame  length  is  the  most  sensitive  to 
buoyancy  effects  from  temperature. 


The  effect  of  buoyancy  on  the  jet  is  determined  by  the  ratio  of  the  buoyancy 
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logjQ(x/d*) 

Fig.  4.12  Scan  of  ATf  for  0  =  7,  i2e  =  20,000,  Uo  =  62.4 m/s,  p  =  2 atm. 


force  on  the  jet  cone  to  the  nozzle  exit  momentum  flux.  Hecalling  the  flame  geometrj’ 
of  Fig.  3.1,  the  buoyancy  force  B  on  the  jet  cein  be  expressed  as 

B  =  g'K'pV  ,  (4.6) 


where  g  is  the  acceleration  of  gra'cit}'.  A/?  is  the  average  change  in  density  owing 
to  heat  release  given  by  Ap  =  —  'p,  and  V'  is  the  volume  of  the  jet  cone. 

The  inertizil  strength  of  a  jet  source  is  normally  characterized  by  the  nozzle  exit 
momentum  flux  Jo ,  given  for  top-hat  exit  velocity  eind  density  profiles  by 


(4.7) 


41 


Becker  &  Yamazziki  (1978)  attempted  to  identify  when  the  jet  is  momentum- 
dominated  and  where  the  transition  to  buoyancy  occurs  by  comparing  the  buoyancy 
force  on  the  jet  cone  to  the  exit  momentum  flux.  The  jet  grows  linearly  with  x  (e.^., 
Eq.  4.21)  so  they  assumed  that  V'  ~  .  They  worked  with  high  heat  release  flames 

and  made  the  further  assumption  that 

Ap  =  Poo  -p  ^  poc  .  (4.8) 

Thus,  they  give  a  characteristic  measure  of  the  buoj'^ancy  force  B  as 

B  =  g  ApV  ^  g  Poo  ■  (4.9) 


Their  resulting  ratio  of  buoyancy  to  momentum  force,  typically  called  the  Richard¬ 
son  number  Ri  is  given  by 


(4.10) 


They  work  with  the  cube  root  of  the  Richardson  number  by  defining  a  nondimen- 
sional  streamwise  coordinate  ^  as 


e  = 


(4.11) 


where 


r,  _  -^(JPood*^  gd* 

Ru  =  - T 

4Jo 


(4.12) 


is  the  source  flame-zone  Richardson  number.  The  source  velocity  Us  is  given  by 


j-  _ 

mo 


(4.13) 


From  Eq.  4.11,  ^  0  corresponds  to  the  momentum-dominated  regime  and  ^  >  oo 

to  the  buoyancy-dominated  regime.  Becker  k,  Yamtizaki  used  a  pitot  probe  and  an 
uncorrected  thermocouple  to  infer  point  measurements  of  density  and  velocity.  The 
mass  flux  in  their  measurements  was  evaluated  numerically  from  the  product  of 
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logjQ(x/d*) 

Fig.  4.13  Normalized  product  thickness  versus  Iogio(x/d*)  for  severed  adiabatic 
flame  temperatures,  Re  =  40,000,  ^  =  10,  J/q  =  62.4 m/s,  p  = 
4  atm. 


these  measurements.  Although  they  express  some  concern  over  the  errors  inherent 

with  their  teclmique,  they  find  that  for  ^  <  1 ,  the  mass  flux  of  the  jet  is  linear  with 

X ,  the  result  expected  from  the  similzLrity  analysis  for  a  momentum-dominated  jet 

(e.^.,  Ricou  Spalding  1961).  For  ^  >  2.5,  their  measurements  indicated  that  the  ® 

mass  flux  vziries  as  ,  which  agrees  with  the  similarity  result  for  a  buoyant  jet 

with  Ap  «  poo  (c.p.,  Delichatsios  ir  Orloff  1984).  Toner  (1987)  performed  more 

direct  entrainment  measurements  utilizing  gas  chromatography  and  found  similar  # 

behavior. 


logjo^  fip/l-w 
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Fig.  4.14  The  data  of  Fig.  4.13  in  log-log  coordinates.  The  -2/3  slope  is  included 
for  reference. 


Becker  &  Yanicizaki  and  Toner  used  high  heat  release  flames  so  their  assumption 
of  Ap  w  Poo  is  probably  valid.  However,  for  the  present  experiments,  Ap  •<  poo 
so  this  assumption  is  not  valid  and  Ap  is  retained  in  the  buoyancy  expression  B 
of  Eq.  4.9.  Assuming  poo  ^  Po  •  '^oc  ^  To ,  that  the  temperature  of  the  entire  jet 
cone  is  scaled  by  Too  +  ATf  (a  conserv’ative  estimate),  and  ideal  geis  behavior  with 
constant  pressure  yields  the  version  of  Eq.  4.11  used  in  the  current  analysis,  :.e.. 
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(4.14a) 


The  Msumption  of  constant  pressure  in  the  jet  cone  is  not  strictly  true.  Townsend 
(1976),  for  example,  notes  that  p  +  v'^  ~  Poo  where  p  is  the  jet  pressure,  v'  is 
the  instantaneous  transverse  velocity  fluctuation,  and  po^  is  the  reservoir  pressure. 
However,  typically  p  for  these  experiments,  so  the  assumption  of  constant 

pressure  is  certainly  valid. 


Experiments  to  test  the  buoyancy  scaling  ^  on  such  parameters  as  flame  tern-  • 

perature  and  exit  velocity  were  performed.  Figure  4.13  shows  a  set  of  experiments 
spanning  adiabatic  flame  temperatures  AT(  from  3.5  K  to  215  K.  A  Reynolds  num¬ 
ber  of  40,000  was  chosen  for  this  sczui.  A  ^  of  10  placed  the  flame  tip  in  the  center  ^ 

of  the  logio(a;/do)  range,  thus  giving  an  equal  number  of  points  on  each  side  of  the 
knee.  The  temperature  measurements  beyond  the  flame  tip  reach  a  constant  value 
for  the  low  heat  release  runs,  indicating  that  they  are  momentiun-dominated.  As  ^ 

the  flame  temperature  increases,  however,  the  temperature  measurements  no  longer 
attain  a  constant  level  but  begin  to  decre2ise  with  increrising  xfd*  .  This  is  evidence 
of  the  jet  begiiming  to  feel  the  effects  of  buoyancy. 


To  compare  the  temperature  decay  after  the  flame  tip  in  Fig.  4.13  with  the 
scaling  discussed  above,  the  data  were  plotted  in  log-log  coordinates  in 
Fig.  4.14.  In  these  coordinates,  a  power  law  behavior  shows  up  as  a  linear  region  # 

with  the  power  of  x  given  by  the  slope.  The  -2/3  slope  is  included  for  reference. 
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Fig.  4.15  Plot  used  to  determine  xt  for  the  ATf  =  109.4 K  run  of  Fig.  4.13. 
Arrow  indicates  the  estimated  location  of  it . 


The  measured  temperatures  after  the  flame  tip  decrease  at  an  increasing  rate  with 
increasing  AT{ .  The  variation  in  the  slope  of  the  data  in  Fig.  4.14  beyond  the  fleone 
tip  from  the  lowest  ATf  to  the  highest  is  roughly  0  to  —1/3.  This  indicates  that  the 
jet  is  in  transition  between  the  momentum-dominated  and  the  buoyancy-dominated 
regime  with  incre2ising  ATf.  This  allows  us  to  conclude  that  the  experiments  of 
Sec.  4.1  were  free  of  buoyemcy  effects;  the  measurements  are  free  of  this  temperature 
decay. 


The  transition  to  the  buoyancy-dominated  regime  in  Fig.  4.13  was  identified 


log  10 

Fig.  4.16  Plot  used  to  detennine  .Tt  for  the  AT{  =  27.6  K  run  of  Fig.  4.13.  Arrow 
indicates  the  estimated  location  of  . 


as  the  point  where  the  temperature  mezisurements  departed  from  the  lowest  ATf 
runs  that  were  momentum-dominated.  By  plotting  the  curve  fits  of  the  two  runs 
with  ATf  =  3.5  and  6.9  K  with  the  data  of  the  other  ATf  nms  individually, 
the  transition  location  Xt  for  each  ATf  was  determined.  Figures  4.15  and  4.16 
demonstrate  the  method  with  the  transition  ^'alue  indicated.  The  error  of  this 
determination  of  it  was  estimated  by  considering  the  location  of  the  closest  wire 
2ind  the  character  of  the  departiue.  Some  of  the  runs  showed  a  more  distinct 
departure  from  the  reference  curves  than  others.  The  reader  is  invited  to  sight 
along  these  curves  to  assess  the  choices  for  it . 
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Fig.  4.17  Values  of  Xti  determined  from  Fig.  4.13,  normalized  by  for  several 
adiabatic  flame  temperatures,  <f>  =  10,  Uo  =  62.4 m/s,  p  =  4 atm. 
Beuts  represent  the  uncertainty  of  the  Xt  estimate. 


To  determine  the  value  of  ^  that  corresponds  to  Xt,  which  we  will  call 
Eq.  4.14a  will  be  rewritten  as  the  product  of  three  nondimensional  groups,  i.e., 


i  = 


(  ATf 

VATf  +  Too/  \hV) 


(4.14b) 


Therefore,  for  constant  To©,  Uo,  and  do,  (  =  fn(Arf,x).  Note  that  there  is  no 
distinction  in  ^  as  to  whether  Xt  is  in  the  flame  zone,  or  beyond  the  flame  tip.  The 
flame  length  was  incorporated  in  Eq.  4.14b  to  investigate  whether  the  functional 
dependence  of  ^  on  ATf  and  Uo  changes  depending  on  whether  Xi  is  less  thzm  or 
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Fig.  4.18  NormaJized  product  thickness  versus  logio{x/d*)  for  several  exit  veloci¬ 
ties,  <^  =  10,  ATf  =  6.9 K. 


greater  than  L{ .  The  new  formulation  is 

/  AT.  \  / 


f  =  (  V  (± 

^  \ATt  +  To.)  VC'oV  Uf 
Solving  for  x/X-f  gives 


(4.14c) 


(4.15) 


If  {  represents  the  correct  scaling  to  describe  the  transition  from  the  momentum- 
dominated  to  the  buoyancy-dominated  regime,  one  should  be  able  to  vary  AT(, 


Xt/L 
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Fig.  4.19  Values  of  Xi,  determined  from  Fig.  4.18,  normalized  by  Z.f„  for  several 
exit  velocities,  </>  =  10,  ATf  =  6.9 K.  Bars  represent  the  uncertainty 
of  the  Xt  estimate. 


holding  Too,  Uq,  and  do  constant,  and  calculate  a  single  transition  value  of  by 
substituting  the  Xt’s  into  Eq.  4.15. 

The  values  of  Xt  that  were  estimated  from  the  data  plotted  in  Fig.  4.13  were 
normalized  by  ,  where  is  the  value  of  Lf  determined  from  the  curve  fit  in 
Fig.  4.6  for  each  particular  Re  and  <l>.  We  can  plot  Xt/Lf^  versus  (1  +  Too/ATf)^^^ 
and  fit  a  line  though  the  data,  forcing  the  fit  to  go  through  the  origin.  The  slope 
from  the  fit  represents  the  coefficient  of  (1  +  Too/ATf)^^^  in  Eq.  4.15.  Using  the 
known  values  of  Uq,  g,  do,  and  L{^ ,  the  value  of  can  be  computed.  The  Xt  data 
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Fig.  4.20  The  data  of  Fig.  4.17  in  logarithmic  coordinates.  The  1/3  slope  is  in¬ 
cluded  for  reference.  See  Fig.  4.17  for  symbol  legend. 


and  the  curve  fit  are  shown  in  Fig.  4.17.  The  bars  indicate  the  uncertainty  of  the 
It  estimate.  The  curve  fit  suggests  a  value  of  =  0.96  ±  0.06. 

A  similar  scan  of  Uq  was  performed  holding  Too,  ^Tf,  and  do  constant.  These 
data  are  plotted  in  Fig.  4.18.  The  collapse  for  the  runs  at  the  same  velocity  is  good 
in  the  ramp  region,  but  the  data  diverge  after  the  flame  tip. 

The  transition  locations  it  in  Fig.  4.18  were  estimated  as  described  above, 
normalized  by  ,  and  plotted  versus  (Uo^/gdo)  in  Fig.  4.19.  The  curve  fit 
results  in  a  value  of  ^t  =  0.99  ±0.13.  This  vedue  coincides  with  the  value  from 
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Fig.  4.21  The  data  of  Fig.  4.19  in  logarithmic  coordinates.  The  1/3  slope  is  in¬ 
cluded  for  reference.  See  Fig.  4.19  for  symbol  legend. 


the  ATf  seem  and  also  with  the  observed  veilues  reported  by  Becker  Sc  Yamazsdei 
and  Toner,  discussed  above. 

Finally,  the  normalized  xt  s  were  plotted  in  the  appropriate  log-log  coordinates 
in  Figs.  4.20  and  4.21  to  check  the  power  law  scaling  of  Eq.  4.15,  and  to  see  if  there 
is  zuiy  cheinge  in  the  power  law  for  x,  <  L{^ .  The  slope  of  the  curve  fits  for  both  are 
close  to  the  1/3  power  law  of  Eq.  4.15.  The  reference  line  is  the  curve  fit  with  the 
slope  fixed  at  1/3.  In  these  coordinates,  Xi  <  Lt^  corresponds  to  negative  values 
of  logio(a;t/I'fo)- 
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We  appreciate  that  the  ^  -scaling  described  above  cannot  be  expected  to  serve 
as  the  correct  description  for  xt  <  L{ .  These  data  were  presented  here  as  a  first  ef¬ 
fort  to  ascertain  that  the  Reynolds  number  study  W2is  conducted  in  the  momentum- 
dominated  flow  regime,  to  provide  a  direct  eissessment  of  buoyancy  effects,  and  to 
investigate  the  role  of  the  various  nondimensional  groups.  The  issue  of  buoyancy 
in  turbulent  jet  flows  is  not  closed,  by  any  means,  and  more  experiments  designed 
for  this  purpose  are  needed  for  a  complete  documentation  of  this  important  phe¬ 
nomenon. 


4.3  Kinetics 

As  noted  earlier,  the  F2/NO  cheniistry  in  this  study  can  be  described  by  the 
“effective”  reaction: 

k-ral 

2N0  +  F2^2N0F,  AQ  =  - 75  --r--  .  (4.16) 

mole  F  2 

The  reaction  actually  consists  of  two  chain  reactions  of  differing  rates  and  heats  of 
reaction  given  by 


NO -H  F2  ^  NOF -t- F,  AQ  =  -18-^^, 

mole 


(4.17a) 


with 

ki  w  4.2x10"  exp^  , 

\  Ka-T  /  mole  •  sec 

where  is  the  universal  gas  constant  and  T  is  the  temperature  of  the  reaction. 
The  second  reaction  is  given  bj^ 

kcal 


F  -h  NO  -f  M  ^  NOF  -j-  M ,  AQ  =  -  57 


mole 


(4.17b) 


with 


3x10'® 


cm 


mole'' 


sec 


AT/AT 
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Fig.  4.22  ChemicaJ  reaction  times  for  some  of  the  rims  of  Sec.  4.1,  p  =  latm, 
T  =  293  K.  See  Table  4.1  for  corresponding  mass  fractions. 


I 


There  is  edso  one  recombination  reaction  given  by 


F  +  F  +  M  ^  F.  +  M , 


AQ  =  -38 


kcal 

mole 


(4.18) 


with 


3.2x10 


14 


cm 


mole'^ 


sec 


I 


The  third  body  M  in  the  reactions  denotes  any  chemical  species  present  in  the 
reacting  mixture.  The  reaction  rates  and  the  heats  of  reaction  AQ  were  taken 
from  Baulch  et  al.  (1981). 
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To  estimate  the  chemical  reaction  time  Tc,  these  reactions  were  investigated 
using  the  CHEMKIN  chemical  kinetics  software  package  (Kee  et  al.  1980)  with  a 
constzint  mziss,  constant  pressure  reactor.  Figure  4.22  shows  a  sample  result  from 
the  progreim  for  the  mass  fractions  corresponding  to  the  four  (j>'s  used  in  Sec.  4.1. 
The  reaction  time  was  determined  by  drawing  a  tangent  to  the  curve  at  its 
steepest  point  and  using  the  intersection  of  this  tangent  line  with  the  adiabatic 
flame  temperature  as  the  chemical  time.  The  slowest  rate  is  for  =  7 ,  p  =  1  atm 
with  Tc  =  0.91msec. 

As  mentioned  earlier,  it  was  required  that  the  time  lag  between  molecular 
mixing  and  chemical  product  formation  be  negligible.  For  this  mixing-limited  con¬ 
dition,  which  we  will  designate  as  the  fast  chemistrj'  regime,  there  is  a  one-to-one 
correspondence  between  molecular  mixing  and  product  formation.  The  ratio  of  the 
local  large-scale  pzissage  time  to  the  chemical  reaction  time  is  one  way  of  estimating 
when  the  reacting  flow  is  in  the  fast  chemistry  regime.  This  ratio  is  the  Dahmkohler 
number  given  by 


Da  = 


n 


(4.19) 


where  rg  is  the  local  large-scale  passage  time,  estimated  by  dividing  the  local  jet 
diameter  S{x)  by  the  mean  centerline  velocity  Ud 


n  = 


Ucl 


(4.20) 


and  Tc  is  the  chemical  reaction  time,  computed  as  described  above.  The  local  jet 
diameter  6(x)  for  a  turbulent  jet  has  been  measured  by  many  experimenters  {e.g., 
White  1974;  Dowling  1988)  and  can  be  approximated  by 


<*>(r)  %  0.44  r  . 


(4.21) 


An  estimate  of  Uc\  is  given  for  a  momentum-dominated,  density-matched  jet  by 
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Chen  &:  Rodi  (1980)  as 

U,x  =  6.2  (4-22) 

The  jet  exit  diameter  do  in  the  above  formula  should  actually  be  replaced  by  d* , 
the  momentum  dieimeter  of  the  jet,  but  for  a  properly  designed  nozzle  with  density- 
matched  jet  and  reser\'oir  gases,  d*  rs  do  (see  Sec.  2.2).  In  the  limit  of  fast  chem¬ 
istry  {Da  1),  chemical  product  formation  is  equivalent  to  molecular  mixing. 
The  lowest  Da  of  the  runs  in  Sec.  4.1  occurs  at  x/d*  =  30  for  the  4>  =  1  run  of 
Fig.  4.22,  with  a  corresponding  Dahmkohler  number  of  Da  =  2.8. 

Since  Uq  is  constant  for  the  Reynolds  number  experiments  described  in  Sec.  4.1, 
and  the  spreading  rate  of  the  jet  is  assumed  to  be  constant,  ts  for  a  given  xjd*  is  the 
same  for  all  of  the  current  experiments  in  Sec.  4.1.  The  chemiceil  time  decreases 
with  reactant  number  density  (i.c.,  pressure).  This  means  that,  for  the  experiments 
described  in  Sec.  4.1,  the  numerator  in  Eq.  4.19  remeuns  constant  and  the  denomina¬ 
tor  is  decreasing  with  increasing  Reynolds  number.  Therefore,  demonstrating  fast 
chemistry  for  the  lowest  Reynolds  number  runs  guarantees  that  the  higher  Reynolds 
number  experiments  of  Sec.  4.1  are  in  the  fast  chemistry  regime. 

One  way  of  determining  if  a  run  is  in  the  fast  chemistry  regime  is  to  repeat 
the  experiment,  doubling  the  concentrations  of  both  jet  and  reservoir  fluids.  This 
doubles  the  adiabatic  flzune  temperature.  The  kinetic  rate  is  proportional  to  reac¬ 
tant  nrnnber  density  and  also  exponentially  dependent  on  the  reaction  temperature. 
Therefore,  a  concentration  doubling  increases  the  kinetic  rate  through  higher  reac¬ 
tant  concentrations  and  by  increasing  the  reaction  temperature.  The  temperature 
rise  for  the  higher  concentration  run  should  double,  if  the  chemical  product  forma¬ 
tion  is  already  mixing-limited.  If  the  chemistry  was  not  fast  before  the  concentration 
doubling,  the  faster  kinetic  rate  after  the  doubling  should  result  in  more  them  dou¬ 
ble  the  chemical  prod’i:t  and,  correspondingly,  more  than  double  the  temperature 


Fig.  4.23  Temperature  doubling  to  check  kinetic  rate,  <j>  =  7 ,  Uq  =  62.4 m/s. 
Dahmkohler  number  is  given  for  x/d*  =  30. 


rise.  Note  that  the  Dahmkohior  numbers  at  x/d*  =  30  in  Fig.  4.13  are  in  the 
range  of  7.8  <  Da  <  1,100,  yet  the  collapse  of  the  data  within  the  ramp  ’•egion 
is  excellent,  indicating  that  Da  >  7.8  is  in  the  feist  chemistry  regime. 

A  simileir  temperature  doubling  was  performed  for  the  run  in  Sec.  4.1  with  the 
lowest  Dahmkohler  number,  shown  in  Fig.  4.23.  These  data  do  not  collapse.  Slow 
chemistry  should  only  affect  the  ramp  region  of  the  curve  and  possibly  result  in 
a  longer  flame  length.  The  asymptotic  level,  however,  should  eventually  agree  at 
some  downstream  location  because  all  of  the  F2  should  react  if  one  goes  far  enough 
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downstream.  The  asymptotic  level  does  not  agree  for  the  various  ATf’s.  This 
raised  the  suspicion  that  the  reactant  concentrations  were  somehow  eiffected. 

Conceivable  inaccuracies  in  the  loading  procedures  could  not  account  for  such 
a  large  drop  in  the  asymptotic  level.  Some  other  mechanism  for  affecting  the  con¬ 
centrations  had  to  be  responsible.  The  fact  that  the  asymptotic  level  nearly  doubles 
for  each  temperature  doubling  lead  to  the  hyjjothesis  that  some  of  the  F2  was  be¬ 
ing  depleted  by  passivation  of  the  nozzle  plenum  screens  and  honeycomb  during 
each  run.  Passivation  is  a  term  commonlj'  used  in  fluorine  systems.  It  refers  to 
the  process  of  admitting  low  concentrations  of  F2  into  the  previously  unexposed 
metal  surfaces  in  a  new  flow  system  to  allow  oxidation  of  the  metal  by  the  fluorine. 
The  oxidation  layer  provides  a  protective  coating  on  the  metal  surface  which  resists 
further  corrosion  by  subsequent  exposure  of  fluorine.  The  <f>  =  7,  ATf  =  7.2  K 
run  at  1  atm  has  the  lowest  concentrations  of  F2  and  NO  of  all  the  runs  in  Sec.  4.1. 
This,  along  with  the  fact  that  this  run  is  also  at  the  lowest  pressure,  results  in  the 
lowest  reactant  partial  pressures  and  therefore  the  least  absolute  number  of  moles 
of  each  reactant.  The  nozzle  is  inside  the  HPRV  and  therefore  sees  NO  during  the 
loading  process  before  the  run.  The  NO  reacts  with  the  F2  passivation  layer  in  the 
plenum  requiring  re-passivation  from  the  F2  in  the  jet  gas  during  the  run.  Assum¬ 
ing  the  depletion  rate  is  constant,  the  effect  would  be  most  pronounced  when  the 
absolute  molm  amount  of  F2  is  lowest,  as  is  the  case  for  the  runs  in  question.  At 
low  concentrations,  doubling  the  absolute  number  of  moles  of  F2  by  doubling  the 
concentrations  should  cut  the  passivation  effect  in  half.  This  appears  to  be  what  is 
happening  in  Fig.  4.23,  supporting  the  passivation  hypothesis. 

The  kinetics  and  the  passivation  problems  disappear  for  the  run  with  ATf  = 
14.4  K  in  Fig.  4.23,  as  demonstrated  by  the  collapse  of  the  data  with  the  higher 
Reynolds  number  run  that  should  be  in  the  fast  chemistry  regime.  For  this  reason, 
the  run  in  Fig.  4.1  at  1  atm  was  done  at  double  the  adiabatic  flame  temperature. 


Fig.  4.24  Temperature  doubling  to  check  kinetic  rate.  Dalunkohler  number  is  given 
for  x/d*  =  30. 


The  collapse  of  the  data  in  Figs.  4.1,  4.2,  4.3,  and  4.4  indicate  that  all  the  runs  of 
Sec.  4.1  are  in  the  fast  chemistry  regime  and  free  of  passivation  problems.  Figure 
4.12  also  proves  that  the  run  at  the  higher  ATf  =  14.4  K  was  still  momentum- 
dominated  because  the  temperature  ine2isurement  beyond  the  flame  tip  remains 
nearly  constant. 

Figure  4.24  shows  the  result  of  a  similar  experiment  at  p  =  2  atm  with  two 
temperature  doublings.  The  collapse  is  good  and  this  condition  appeeirs  to  be  in 
the  fast  chemistry  regime  with  the  A  Tf  =  7.2  K  run  having  a  Da  =  5.6 . 
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logjo(x/d*) 


Fig.  4.25  Runs  showing  the  lowest  Dahmkohler  number  that  still  demonstrates  feist 
chemistry.  Da  is  given  for  x/d*  =  30. 


Figure  4.25  compares  the  run  at  1  atm,  4>  =  14,  to  the  run  at  15  atm  and 
the  same  (j)  and  flame  temperature.  This  is  the  lowest  value  of  the  Dalimkohler 
number  in  all  of  the  experiments  which  still  collapses  well  with  an  experiment  at  a 
much  leu-ger  Dahmkohler  number.  The  lower  limit  of  Dahmkohler  number  in  these 
experiments,  in  which  the  flow  can  be  considered  to  be  in  the  fast  chemistry  regime, 
is  therefore  Da  ~  5.2 . 

In  the  runs  described  thus  far,  the  Dahmkoler  number  was  clianged  by  varying 
the  chemical  time,  but  changing  the  large-scale  passage  time  is  another  way  of 
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Fig.  4.26  Runs  showing  the  effects  of  sIo%v  chemistry  in  the  reimp  region,  =  10, 
ATf  =  6.9  K.  Dahmkohler  number  is  given  for  x/d*  =  30. 


varying  Da.  Several  runs  at  different  velocities  emd  pressures  were  performed  in 
which  the  Dalimkohler  number  was  low  but  the  temperature  measurement  beyond 
the  flame  tip  achieved  the  correct  level,  showing  them  to  be  free  of  passivation 
problems.  Figure  4.26  shows  runs  at  various  velocities  and  Dahmkohler  numbers. 
The  temperatures  for  the  slow  chemistry  conditions  lag  behind  the  fast  chemistry 
runs  but  eventually  achieve  the  same  level  beyond  the  flame  tip.  This  is  how  slow 
chemistry  was  expected  to  manifest  itself.  The  Dalmikohler  numbers  are  consistent 
with  the  limiting  veilue  of  Da  ~  5.2  discussed  above.  The  lower  velocity  runs  show 
the  effects  of  buoyancy  bej’ond  the  flame  tip.  The  high-velocity  rim  is  repeatable. 
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CHAPTER  5 

Conclusions 


A  new  high  pressure  combustion  facility  has  been  built  to  investigate  mix¬ 
ing  in  axisymmetric,  gas-phase  turbulent  jets.  The  study  was  aimed  at  identifying 
Reynolds  number  effects  on  turbulent  mixing  by  examining  changes  in  flame  length, 
which  is  related  to  the  molecular  mixing  rate,  with  increcising  Reynolds  number. 
The  run  parameters  were  chosen  so  that  the  jet  remained  momentum  dominated 
to  the  farthest  mecisuring  station  and  that  the  time  lag  between  molecular  mix¬ 
ing  and  chemical  product  formation  was  negligible.  The  facility  uses  fluorine  and 
nitric  oxide,  diluted  with  nitrogen,  to  provide  a  feist,  hypergolic  chemical  reaction 
with  heat  release.  The  average  temperature  was  recorded  with  a  set  of  long,  thin, 
resistance  wire  thermometers  stretched  across  the  jet  centerline  at  16  downstream 
locations  ranging  from  x/do  =  30  to  240.  Runs  at  several  stoichiometric  mixture 
ratios  were  performed  at  Reynolds  numbers  ranging  from  10, 000  <  Re  <  150, 000 
to  determine  the  flame  length  dependence  on  Reynolds  number.  The  conclusions 
are  summarized  below. 

1.  The  time-averaged  line  integral  of  temperature  performed  by  the  cold  wires 
along  the  transverse  axis  of  the  jet  displays  a  logarithmic  dependence  on 
x/d*  within  the  flame  zone.  Beyond  the  flame  tip,  this  integral  asymptotes 
to  a  constant  value,  as  predicted  from  the  scaling  and  similarity  laws  for 
a  momentum-dominated  jet. 
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The  following  two  conclusions  represent  the  major  findings  of  this  investigation, 
namely  that  flame  length  Li  is  a  function  of  Reynolds  number. 

2.  Runs  at  several  stoichiometric  mixture  ratios  for  increasing  Reynolds  num¬ 
bers  indicate  that  the  normalized  flame  length  Li/d*  displays  a  linear 
dependence  on  and  that  the  slope  A  of  the  linear  fit  decreeises  from 
Rt  —  10,000  to  20,000.  Dowling  h  Dimotakis  (1990)  observed  a  simi¬ 
lar  decrease  in  their  concentration  measurements  of  the  same  magnitude 
over  a  similar  Reynolds  number  range.  The  slope  then  remains  constant 
with  a  value  of  j4  =  5.6  for  20,000  <  Re  <  150,000.  This  threshold 
Rejmolds  number  of  20.000  is  close  to  the  value  of  Re  =  25, 000  found  by 
Ricou  &  Spalding  (1961),  where  their  entrainment  measurements  became 
independent  of  Reynolds  number. 

3.  The  flame  length  measurements  revealed  the  existence  of  a  “mixing  virtual 
origin”  Lq  ,  defined  as  the  far-fleld  flame  length  extrapolated  to  <i>  =  0 . 
Alternatively  stated,  Lq  is  the  flame  length  that  would  result  in  the  limit 
of  </)  — ♦  0  if  the  flow  could  be  described  by  its  far-field  behavior  at  the 
nozzle  exit.  This  Lq  has  a  nonzero  value  that  increases  up  to  Re  —  20, 000 
and  then  decreases  up  to  the  highest  Reynolds  number  investigated  in  the 
experiments  reported  here,  he..  Re  =  150,000. 

4.  Experiments  to  determine  the  parameter  range  where  the  jet  is  mo¬ 
mentum-dominated  were  performed  and  it  was  demonstrated  that  the  con¬ 
clusions  above  were  based  on  runs  that  were  in  the  momentum-dominated 


regime. 
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5.  Experiments  to  identify  the  location  it  where  the  transition  to  the 
buoyancy-dominated  regime  occurs  indicate  that  the  nondimensional 
downstream  buoyancy  coordinate  ^  works  reasonably  well  in  describing 
this  point  for  a  variety  of  run  conditions,  for  it  >  Li .  The  transition 
value  was  found  to  be  ^t  —  1-  This  value  agrees  with  Becker  k.  Ya- 
mazaki  (1978)  and  Toner  (1987),  who  used  entrainment  measurements  to 
determine  the  transition  location. 

6.  The  chemical  reaction  in  the  jet  was  found  to  be  in  the  fast  chemistry 
regime  for  Da  >  5.2.  All  of  the  runs  on  which  the  above  conclusions  were 
based  were  shown  to  be  in  the  fast  chemistry  regime. 


APPENDIX  A 


Experimental  design  and  run  parameters 

The  High  Pressure  Combustion  Facility  (HPCF)  was  designed  to  perform  jet 
mixing  experiments  over  a  pressure  range  of  0.1  atm  to  15  atm  using  both  conven¬ 
tional  fuels  and  more  exotic  reactants  such  «is  fluorine,  nitric  oxide,  and  hydrogen. 
The  lab  building  is  a  4.57  m  by  5.79  m  (15  ft  by  19  ft)  extension  to  the  HF  Shecir 
Layer  Combustion  Facility.  The  list  of  desired  run  parameters,  shown  in  Table  A.l, 
determined  much  of  the  design  of  the  facility.  The  following  sections  give  brief  de¬ 
scriptions  of  the  major  experimental  items  with  some  of  the  motivating  factors  and 
a  short  description  of  the  run  parameters  and  procedures. 


A.l  High  Pressure  Reactant  Vessel  (HPRV) 

The  tank  specifications,  shown  in  Table  A. 2,  resulted  from  many  constraints. 
The  diameter  was  determined  by  the  requirement  of  low  recirculation  velocity  of  a 
confined  jet  {Ui  <  2%  of  L’cl  at  ^/do  =  240)  and  the  height  was  set  to  be  one  half 
of  the  floor  to  ceiling  height.  The  windows  are  positioned  10.16cm  (4  in.)  above 
the  vessel  centerline  to  maximize  the  jet  axial  station  that  can  be  viewed.  There 
are  44  internal  stainless  steel  tabs  for  use  in  supporting  internal  hardware.  Fluorine 
is  extremely  corrosive  so  the  choice  of  tank  materials  was  either  stainless  steel  or 
mild  steel  with  a  teflon  coating.  Stainless  steel  was  cost  prohibitive  so  the  tank  was 
made  of  mild  steel  with  an  internal  coating  of  127 ptm  (.005  in.)  thick,  pinhole-free 
PFA  teflon,  which  has  a  maximum  operating  temperature  of  644  K.  The  temk  is 
mounted  with  the  base  supported  1.22m  (4  ft)  above  the  floor  and  the  top  0.61m 
(2  ft)  from  the  ceiling  (see  Fig.  2.2).  This  allows  the  lower  head  to  be  dropped  down 
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Table  A.l:  High  Pressure  Combustion  Facility  Target  Specifications 

ITEM  DESCRIPTION 


Reservoir 

•  Reactant  Gases 

•  Diluent  Gases 

•  Temperature  Range  (w/o  liner) 

•  Temperature  Range  (w/  liner) 

•  Concentration  Range  (molar) 

Jet 

•  Reactant  Gases 

•  Diluent  Gases 

•  Temperature  Range  (unreacted) 

•  Nozzle  Exit  Diameter 

•  Maximum  Mass  Flow  Rate 

•  Gas  Delivery  Method 

•  Concentration  Range  (molar) 

Overall 

•  Maximum  Run  Time 

•  Reynolds  Number  Range 

•  Pressure  Range 

•  Flame  Temperature  (Fj/NO) 

•  Flame  Temperature  (Hydrocarbons) 


H2 ,  NO,  air 

N2,  Ar,  He 

300  K  <  T  <  644  K 

300  K  <  T  <  2,500K 

0  <  Chj.no  <  8% 


F2 ,  hydrocarbons 
N2,  Ar,  He 
Ambient  (300K) 

2.5  mm 

0.01  kg/s  (N2) 

Expulsion  Diaphragm 
0  <  Cf,  <  40% 

68  seconds  (At  Uq  =  125 m/s) 
1,000  <  Re  <  150,000  (nominal) 
0.1  atm  <  p  <  15  atm 

3.5  K  <  ATf  <  215  K 
ATf  ~  2500  K 
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Table  A. 2:  High  Pressure  Reactant  Vessel  Specifications 


Manufacturer: 
Serial  No.  : 
Year  Built  : 


California  Tank  and  Manufacturing  Corporation 

21274 

1985 


Working  Pressure  @  Temp:  full  vacuum  to  15  atm  ^  422  K 
Hydrostatic  Test  Pressure:  21.4  atm 


Material  :  SA-455-75  steel  (shell  and  both  caps) 

Thickness  :  9.525mm  (0.375  in.)  (shell  and  caps) 


Overall  Length  : 
Nominal  O.D.  : 

Caps  ; 

Internal  Capacity  : 
Overall  Empty  Weight: 


183  cm  (72  in.)  (does  not  include  welding  neck  flanges) 
107  cm  (42  in.) 

91.44  cm  (36  in.)  crown  radius; 

8.89  cm  (3.5  in.)  knuckle  radius 
1.43  m^  (50.Sft^) 

2,268kg  (5000 lbs)  approx. 


Openings  :  8  ea  -  7.62  cm  (3  in.)  welding  neck  flanges  (cap) 

2  ea  -  15.24  cm  (6  in.)  welding  neck  flanges  (cap) 

8  ea  -  7.62cm  (3  in.)  welding  neck  flanges  (shell) 

4  ea  -  3.81  cm  (1.5  in.)  welding  neck  flanges  (shell) 
1  ea  -  40.64cm  (16  in.)  manhole 

3  ea  -  25.4cm  (10  in.)  diameter  windows  (shell) 


Notes  : 


IntPinal  127 /im  (0.0C5  in.)  PFA,  pinhole-free 
teflon  coating  applied  by  Thermech  Engineering. 
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and  the  upper  head  to  be  lifted  free  of  the  shell  giving  full  interior  access.  With 
this  arrangement,  the  tank  windows  are  located  2.13m  (7  ft)  above  the  floor,  a 
desirable  safety  feature  in  case  of  window  blowout  and  for  future  laser  beam  access. 


A. 2  Fluorine  Reactant  Tank 

The  jet  gas  is  contained  in  a  stainless  steel  reactant  tank  with  an  internal 
volume  of  0.042  m^  ( 1 .48  ft^  ).  A  teflon  expulsion  diaphragm  separates  the  jet  charge 
from  the  inert  driver  gas.  The  reactant  tank  is  internally  lined  with  stainless  steel 
screen,  both  to  support  the  diaphragm  over  the  inlet  and  exit  ports,  and  to  provide 
a  gas  path  to  prevent  trapped  pockets  of  gas  during  evacuation  and  expulsion.  The 
reactant  tank  is  located  underneath  the  HPRV  to  minimize  the  length  of  the  delivery 
line  to  the  nozzle,  resulting  in  the  shortest  possible  startup  times  (~  0.8 sec). 

The  teflon  expulsion  diaphragms  were  fabricated  at  Thermech  Engineering  of 
Anaheim  from  a  highly  polished,  hard-chrome  plated  male  mold  made  at  Caltech’s 
Central  Engineering  Services.  The  process  involves  charging  the  mold  with  a  static 
charge,  spraying  it  with  a  layer  of  PFA  teflon  powder,  and  then  baking  it  in  an  oven 
at  700  K  until  the  powder  melts  and  flows  to  form  a  smooth  layer  approximately 
76/.im  (0.003  in.)  thick.  The  process  is  repeated  until  the  desired  thickness  of 
203 //m  (0.008  ir.)  is  obtained.  The  diaphragm  was  periodically  inspected  after  a 
number  of  runs  and  replaced  when  necessary. 


A. 3  Gas  Delivery  and  Exhaust 

Designing  a  gas  deliver}’  system  to  handle  F2  and  NO  gases  involved  many 
special  design  features.  Fluorine  is  extremely  reactive  while  NO,  in  the  absence  of 
water  vapor  or  oxidizers,  is  fairly  inert.  Both  gases  are  toxic  with  potentially  lethal 
concentrations  as  low  as  1  ppm,  necessitating  a  gas-tight  environment.  The  issue 
of  safety  required  that  all  exposed  surfaces  be  well  cleaned  and  the  gas  paths  kept 
leak-tight  and  free  of  water  vapor.  To  ensure  cleanliness,  the  F2  reactant  tank  was 
nitric  acid  pickled,  the  F2  mixing  vessels  were  electropolished,  and  the  NO  mixing 
vessels  were  clear  anodized  prior  to  installation.  Before  any  runs  were  performed, 
the  gas  lines,  the  HPRV,  the  F2  reactant  tank,  and  the  mixing  vessels  were  cleaned 
with  an  M-17  degreaser,  purged  with  dry  nitrogen,  sealed,  and  evacuated.  The  F2 
lines  and  F2  reactant  tank  were  then  exposed  to  low  concentrations  of  F2  for  a 
period  of  time  (typically  2  hours)  to  build  up  a  fluorine  passivation  layer  on  the 
metal  surfaces.  The  entire  system  was  helium  leak-checked  with  a  Matheson  helium 
leak  detector  ( 10~®  cc/min  resolution).  To  ensure  that  no  atmospheric  water  vapor 
enters,  the  system  is  always  kept  at  some  positive  gauge  pressure  with  dry  nitrogen. 
To  remove  any  contaminants  or  water  vapor  that  could  be  present  in  the  gas  supply, 
all  the  nitrogen  used  in  the  experiment  was  passed  through  Robbins  Aviation  13X 
molecular  sieves. 

The  reactive  nature  of  Ft  limited  the  choice  of  materials  for  the  F2  delivery 
system,  ^’irtually  all  of  the  gas  tubing  and  valves  for  the  F2  gas  system  are  stainless 
steel  because  of  its  high  corrosion  resistance.  The  F2  mixing  vessels  are  31GL 
stainless  steel  with  removable  flanges  sealed  with  lead  gaskets.  Although  NO  is 
not  as  reactive  as  F2  ,  the  gas  tubing  and  valves  for  the  NO  delivery  system  are 
also  stainless  steel.  The  NO  mixing  vessels  are  606TT6  aluminum  with  removable 
flanges  sealed  with  \bton  0-rings. 


The  dynamic  ressure  of  the  nozzle  for  the  entire  range  of  experiments  varied 
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from  2  torr  to  500  torr.  To  maximize  the  resolution  of  this  pressure  measurement, 
three  Barocel  differential  pressure  tremsducers  of  reinges  10,  100,  and  1000  torr 
were  available.  The  dynamic  pressure  of  the  run  determined  which  of  the  three 
transducers  would  be  used  to  obtain  the  best  resolutioii. 

There  are  three  distinct  exhaust  systems  for  the  facility.  The  first  is  a  Camfiex 
relief  valve  installed  on  the  HPRV  to  relieve  any  pressure  rises  that  occur  during 
runs  at  pressures  >  1  atm .  The  exhaust  from  this  valve  is  plumbed  into  the  down¬ 
stream  catch  bag  shared  with  the  HF  Shear  Layer  Combustion  Facility  described 
in  Hall  (1991).  The  bag  material  is  J-22,  a  vinyl-based  polymer  that  was  tested 
for  NaOH  and  F2  compatability  in  the  HF  Shear  Layer  Facility.  The  bag  has  a 
capacity  of  120  m^  (4,23Sft^)  and  is  equipped  with  an  internal  NaOH  and  water 
shower  system  to  neutralize  the  exhaust  products.  The  second  system  is  a  132-liter 
(35-gallon)  stainless  steel  barrel  filled  with  NaOH  and  water.  Any  imreacted  NO 
in  the  mixing  vessels  and  by-products  left  in  the  HPRV  after  a  run  are  bubbled 
through  the  barrel  to  neutralize  these  chemicals.  The  third  is  a  20S-liter  (55-gallon) 
stainless  steel  barrel  filled  with  activated  charcoal  in  series  with  a  132-liter  (35- 
gallon)  stainless  steel  drum  filled  with  NaOH  and  water.  Similarly,  the  unreacted 
F2  in  the  mixing  vessel  and  reactant  tank  is  passed  though  the  charcoal  barrel  and 
bubbled  through  the  NaOH  barrel  to  neutralize  the  gas  (see  Fig.  2.1  for  a  schematic 
of  this  system). 

Relief  of  accidental  overpressure  in  the  HPR\'^  and  the  NO  mixing  vessels  is 
provided  by  Fike  burst  diaphragms  installed  on  each  vessel.  The  exhaust  from 
these  diaphragms  is  plumbed  into  the  downstream  catch  bag.  Burst  diaphragms 
are  also  installed  on  the  F2  mixing  vessels  and  a  relief  valve  is  installed  on  the  F2 
reactant  tank,  all  plumbed  to  the  charcoal/NaOH  exhaust  system  described  above. 
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A. 4  Run  Parameters 

The  design  of  the  HPCF  was  based  on  a  nominal  run  condition  of  N2  at 
latm,  with  Uq  =  60  m/s  and  do  =  0.005  m,  resulting  in  a  nominal  Re  of  20,000. 
The  various  components  were  then  sized  to  proi'ide  run  times  sufficiently  long  to 
extract  reliable  statistics,  i.e.,  before  recirculation  of  hot  products  contaminated  the 
temperature  measurements.  To  determine  the  statistical  reliability  of  the  average 
temperature  measurement,  a  local  large-scale  passage  time  can  be  estimated 
from  Eq.  4.20.  This  estimate  can  be  used  to  determine  the  number  of  local  large- 
scale  passage  times  sampled  during  a  run  at  a  particular  x/do  location.  At  the 
nominal  run  conditions  given  above,  for  an  estimated  run  time  of  10  seconds,  this 
corresponds  to  59  's  at  x /do  =  240 . 


The  local  large-scale  passage  time  can  be  compared  to  the  estimated  response 
time  of  the  cold  wires  to  determine  the  minimum  temporal  resolution  require¬ 
ments.  The  response  time  of  the  wires  can  be  estimated  from  an  energy  balance 
between  an  elemental  length  of  wire  and  the  gas  flowing  over  the  wire,  neglecting 
conduction,  radiation,  and  ohmic  heating  of  the  wire.  The  resulting  response  time 
is  given  by 

dvi  Pw  Cw 


4  Nu  ke 


where 


dw  =  cold  wire  diameter 
Pv,-  =  cold  wire  density 
c,v  =  specific  heat  of  the  cold  wire 
Nu  =  Nusselt  number 
^’g  =  thermal  conductivity  of  gas. 

Since  Nu  ~  y/We  {e.g.,  Kramers  1946),  the  response  time  of  the  wire  varies  in¬ 
versely  with  \/Rt. 
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Table  A. 3:  Time  scales  for  Re  =  10,000  experiments, 

Uq  =  62.4 m/s,  do  =  2.5mm,  ATf  =  6.9  K 

x/do 

Uc\  (m/s) 

Tf,  (msec) 

Tw  (msec) 

#  ’s  sampled 

29.8 

12.98 

2.52 

5.83 

1,984 

34.0 

11.38 

3.29 

6.16 

1,520 

39.4 

9.81 

4.41 

6.56 

1,134 

46.2 

8.37 

6.07 

7.01 

823 

52.4 

7.38 

7.81 

7.39 

640 

68.4 

5.66 

13.30 

8.23 

376 

80.2 

4.82 

18.30 

8.76 

273 

91.0 

4.25 

23.55 

9.20 

212 

104.8 

3.69 

31.23 

9.72 

160 

120.0 

3.22 

40.94 

10.23 

122 

138.6 

2.79 

54.62 

10.79 

91 

159.0 

2.43 

71.88 

11.35 

70 

183.2 

2.11 

95.43 

11.94 

52 

208.4 

1.86 

123.48 

12.49 

40 

240.2 

1.61 

164.04 

13.11 

30 

After  performing  runs  with  a  5  mm  exit  diameter  nozzle,  it  was  apparent  that 
the  recirculation  problem  was  worse  than  estimated  and  the  only  solution  was  to 
enlarge  the  HPRV  or  provide  coflow  to  satisfy  the  entrainment  requirements  of  the 
jet.  Since  neither  of  these  could  be  simply  implemented,  the  nozzle  exit  diameter 
was  cut  in  half  to  essentially  “shrink”  the  flow  and  make  the  HPRV  effectively 
larger  in  the  sense  of  increasing  uncontaminated  run  times.  This  set  the  nominal 
run  condition  to  N2  at  1  atm  with  Uq  =  62.4 m/s  and  do  =  2.5mm,  resulting  in 
a  nominal  Re  of  10,000.  With  the  pressure  range  of  0.1  atm  to  15  atm,  this  gives  a 
nominal  Reynolds  number  range  of  1000  <  Re  <  150,000.  The  actual  Reynolds 
numbers  realized  for  the  current  experiments  (at  the  nominal  exit  velocity  and 
diameter)  were  10,000  <  Re  <  150,000.  The  lower  limit  was  set  by  the  F2 
passivation  problem  discussed  in  Sec.  4.3,  while  the  upper  limit  was  set  by  the 
pressure  limitations  of  the  facility. 


The  relevant  times  for  the  lowest  Reynolds  number  runs  are  given  in  Table  A. 3. 
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The  wire  response  times  in  Table  A. 3  are  conservative  because,  as  shown  above,  the 
frequency  response  of  the  wires  increases  with  increasing  Re.  For  this  Reynolds 
number,  the  measurements  should  become  time- resolved  for  x/do  >  52.4.  However, 
since  only  mean  quantities  are  required,  the  temporal  resolution  is  adequate  when 
compared  to  the  number  of  ’s  sampled  during  a  run  to  give  an  accurate  average 
temperature  measurement  at  all  the  x/do  positions.  For  the  runs  where  Uq  was 
varied,  similar  calculations  were  made  to  confirm  the  statistical  accuracy  of  the 
mean  temperature  measurement. 


A. 5  Run  Procedures 

Because  of  the  hazardous  nature  of  the  gases  used  in  the  experiment,  a  lengthy 
set  of  written  procedures  are  followed  for  each  run,  with  two  people  always  present. 
A  typical  run  and  cleanup  operation  requires  most  of  a  day.  The  run  procedure  is 
as  follows; 

1.  Wire  calibration 

2.  Evacuation  of  the  HPRV,  F2  reactant  tank,  and  mixing  vessels 

3.  Loading  of  NO  and  N2  into  the  HPRV  to  the  designated  run  pressure 

4.  Loading  of  F2  and  Nt  into  the  F2  reactant  tank 

5.  Wait  a  period  for  the  gases  to  homogenize 

6.  Perform  the  rim  with  data  acquisition 

7.  Unload  the  catch  bag,  HPRV,  F2  reactant  tank,  and  mixing  'essels. 
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APPENDIX  B 

Calibration  and  processing  techniques 


B.l  Cold  wire  calibration 

In  order  to  convert  the  voltage  output  of  the  cold  wires  into  a  temperature 
measurement,  calibration  is  needed.  One  method  is  to  immerse  the  wires  in  a 
gas  of  known  temperature  until  equilibrium  with  the  gas  is  achieved,  record  the 
voltage,  then  change  the  temperature  of  the  gas  and  repeat  the  measurement.  The 
conversion  constants  can  then  be  computed  assuming  a  linear  relation  between 
voltage  and  temperature. 

The  temperature  change  for  the  Cedibration  is  accomplished  by  an  isentropic 
expansion  of  the  gets  inside  the  HPRV.  The  method  is  to  overpressure  the  HPRV  the 
night  before  a  run,  and  let  the  gas  come  to  thermal  equilibrium  overnight.  Great 
care  is  taken  to  eliminate  thermal  gradients  in  the  lab  by  running  recirculating 
fans  continuously.  In  the  morning,  the  wire  voltages  are  recorded  on  the  computer 
during  the  entire  calibration  process  that  lasts  about  seven  seconds.  This  involves 
beginning  the  data  acquisition  about  one  second  before  the  expansion  to  obtain  the 
voltage  offsets,  opening  the  Camflex  valve  (typically  one  second),  and  waiting  about 
five  more  seconds  for  the  HPR\’  to  exhaust  the  gas  to  atmospheric  pressure.  In  this 
way,  all  IS  cold  wires  are  simultaneously  exposed  to  a  uniform  temperature  drop. 
The  temperature  drop  is  computed  from  the  HPRV  pressure  trace  assuming  an 
isentropic  expansion  of  the  gas  and  the  calibration  constants  are  determined.  Figure 
B.l  shows  a  sample  calibration  of  the  18  cold  wires.  The  calibration  constants  do 
not  vary  by  more  than  1%  from  day  to  day  and  are  fairly  constant  even  after  many 


t  (sec) 

Each  channel  offset  by  0.25 K 


Fig.  B.l  Sample  cold  wire  traces  from  a  calibration  temperature  drop.  The  lower 
trace  is  the  reference  cold  wire  at  i/do  =  0,  followed  by  the  16  cold 
wires  from  x/dg  =  30  to  240,  and  finally,  the  other  reference  cold  wire 
at  I /do  =  240. 


hot  runs.  Since  the  ccJibration  is  performed  using  a  temperature  drop  and  the  wires 
measure  temperature  rises  during  a  run,  it  was  necessary  to  check  for  hysteresis.  The 
temperature  rise  in  the  HPRV  during  a  run  is  computed  from  the  HPRV  pressure 
trace  assuming  an  isentropic  compression  and  compared  to  the  temperature  rise 
recorded  by  the  reference  wire  at  x/do  =  0.  The  two  measurements  always  show 
excellent  agreement. 
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B.2  Processing  techniques 

Twenty  channels  of  data  were  recorded  during  a  run:  the  16  cold  wires,  the 
two  reference  cold  wires,  the  HPRV  pressure,  and  the  nozzle  dynamic  pressure. 

After  a  run,  the  data  are  processed  as  follows.  The  data  were  transferred  to 
the  HYDRA  VMS/ VAX  cluster  for  processing.  A  program  calciilated  the  jet  and 
reservoir  gas  mixture  properties  based  on  Thompson  (1984),  and  then  passed  a 
smoothing  filter  over  the  raw  data  to  compress  the  record.  The  offset  voltage  of 
each  channel  was  removed  and  the  wire  voltages  were  converted  to  temperature  rises 
using  the  calibration  constants  just  described.  Using  the  gas  mixture  properties, 
the  nozzle  differential  pressure  was  converted  to  velocity  and  the  HPRV  pressure 
rise  was  converted  to  a  temperature  rise  assuming  an  isentropic  compression. 

The  temperature  rise  from  the  16  cold  wires  and  two  reference  wires  along  with 
a  scaled  velocity  trace  were  then  plotted  as  in  Fig.  3.2.  A  suitable  averaging  region 
was  chosen,  beginning  at  the  time  where  the  velocity  had  reached  steady  state  and 
ending  at  the  time  when  the  recirculated  hot  products  reached  the  reference  wire 
at  x/do  =  240.  The  average  temperature  rise  of  each  wire  was  then  computed  emd 
the  average  temperature  in  the  HPRV,  determined  from  the  HPRV  pressure  trace, 
was  subtracted  from  the  wire  temperatures.  Then,  the  average  wire  temperatures 
were  normalized  by  the  adiabatic  flame  temperature  rise.  The  nondimensional 
product  thickness  6p/L^,  defined  by  Eq.  3.5,  along  with  the  corresponding  x/do, 
was  written  to  a  data  file  where  the  momentum  diameter  was  implemented.  The 
data  were  then  plotted,  as  shown  in  Fig.  B.2,  for  example.  A  least-squares  fit  was 
done  in  the  ramp  region  and  a  cubic  spline  was  fit  to  the  knee  region  (the  solid 
curve).  The  first  derivative  of  the  curve  fit  was  then  plotted  and  used  to  locate  the 
maximum  Sp  /Lw  •  Finally,  a  99%  of  maximum  6p  / Lw  was  determined  and  used  to 
define  the  flame  lengths  discussed  in  Ch.  4. 
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